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I. PREFACE 

A thorough investigation encompassing theoretical study, experi- 
ment and systems analysis and founded upon an existing treatment of 
the oxidation of materials was made of ATJ graphite, a well charac- 
terized basic material, and a Material 11 XV , a reactor material of 
special interest. The existing continuum theory is extended to cover 
non -continuum situations for graphitic materials. Based on approx- 
imate analytic relations, oxidation behavior is predicted for both 
reaction rate controlled and diffusion controlled oxidation for low 
Reynolds number continuum flow and molecular flow. Relations in 
the transition region between these two regimes are also determined. 
Behavior of the reacting surface is described in terms of drag, heat 
transfer and mass transfer for these flow conditions, which corres- 
pond to very high altitudes. The results of the analytical investigation 
are presented in forms suited to engineering calculations. They 
demonstrate that as a non-continuum environment develops, oxidation 
behavior departs from that of continuum. Further, the differences 
that are predicted depend upon whether oxidation is reaction rate or 
diffusion controlled and can be correlated with the shock Reynolds 
number corresponding to the flow conditions existing in the environment. 
The necessary correlating expressions are established and are pre- 
sented in graphic form. 

The experimental investigation provided extension and corroboration 
of existing experimental results based on the continuum theory. 
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Additional experimental data were acquired through a range of condi- 
tions other than high Reynolds number continuum. The results of the 
experiments substantiated the theoretical findings and provided insight 
into the shape change behavior of spherical specimens of ATJ graphite 
and the designated Material "X". Extensive oxidation mass loss data 
were acquired for both materials. Values of reaction coefficients and 
activation energies were determined from the experimental data and a 
positive corroboration of the diffusion controlled regime was established 
for ATJ graphite. In the case of Material n X', the experimental results 
did not clearly define a region of constant mass rate corresponding to 
diffusion controlled oxidation. An estimated value for the diffusion con- 
trolled oxidation of this material was not confirmed. 

Shape changes observed on spherical specimens were, in general, as 
anticipated. Locally, mass loss was dependent upon the local heat trans- 
fer rate and spin. In effect, spin reduces the local rates of heat transfer, 
resulting in lower mass loss on rotating bodies. In the case of Material 
11 X", some peculiar transient variations in mass loss behavior were ob- 
served, as compared with graphite. 

Using the analytical tools and appropriate experimental data, the re- 
entry of various sized spheres and axially oriented cylinders of Material 
"X" on specific trajectories was investigated by programmed machine 
computation. The basic trajectories, treated for bodies with variable 
W/Cj-yA, had parameters into which perturbations were introduced to study 
performance changes. The computations determined that all of the 
particle sizes investigated would readily survive the selected re-entry 
environments. 
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Computations of oxidation performance under conditions approxi- 
mating the test facility environments were carried out for several 
specimens. The results offer a comparison of theoretically founded 
computations with direct experimental evidence. It was found that 
mass losses and dimensional changes showed quantitatively mixed 
agreement, however the trends are quite well defined. Introduction of 
a suitable spin factor into the machine computation to predict losses 
from rotating specimens was necessary to provide comparison with 
the experimental results. The determined spin factor was satisfactory 
in that the computer data showed discrepancies completely analogous to 
those of the fixed sphere, when compared to the experimental data. Since 
the discrepancies can be adequately accounted for upon realization 
that the flight theory equations are only generally applicable to the 
description of the flow properties produced in an arc heated facility, 
it is fairly obvious that further iteration of the computer analyses to 
account for these differences could with relative ease converge the 
disparities. This was accomplished for the case of the stable axially 
oriented cylinder. 

Additionally, some empirical correlations were determined from the 
results of tests on fixed and tumbling short cylinders such that useful 
comparisons can be made with respect to particle shape and spin rate. 
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II. THEORETICAL CONSIDERATIONS 

A. INTRODUCTION 

In describing the entry burning of an object into the earth's 
atmosphere, it is necessary to describe three distinct flight regimes. 

At altitudes above 450, 000 feet the object can be considered to be in the 
free molecule flow regime; a regime in which one can ignore inter - 
molecular collisions. As the object descends towards lower altitudes, 
particles which have been reflected from the vehicle will begin to collide 
with the hyperthermal free stream particles. The effect of these inter- 
molecular collisions will be to steepen density gradients and so give rise 
to the formation of a shock wave. This is referred to as the transition 
regime since it delineates the onset of the region in which macroscopic 
gas phenomena begin to dominate over microscopic phenomena in the 
bulk gas flow. At still lower altitudes, the intermolecular collisions 
become dominant and the physico-chemical phenomena that are present 
are best described by the equations of continuum gas dynamics -- this 
is referred to as the continuum regime . 

An objective of this report is to give approximate 
analytic relations for the prediction of drag, heat transfer, and mass 
transfer over the entire spectrum of entry flight regimes for an ablating 
graphite object. These relations are presented as correlation equations 
which span the free molecular to continuum regimes, and are in a form 
which are quite accessible to engineering calculation. In the following 
sections, the formulation of these correlation relations are discussed, 
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the approximations and assumptions are stated and the equations are 
presented in a form which facilitates their use in engineering calcula- 
tions. 
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B. ANALYSIS OF AEROTHERMOCHEMICAL INTERACTIONS IN 
LOW DENSITY HYPERSONIC FLOW 

1* Approach 

In predicting the entry behavior of small graphite particles it is 
of critical importance to understand the oxidation processes that occur 
under conditions of rarefied gas flow. This is because the smaller 
particles will spend a significant portion of their entry time at the 
higher altitudes where the mean free path for intermolecular collisions 
will be large comparable to the particle size. T o assess the effects of oxida- 
tion in these rarefied flow regimes it is, therefore, necessary to analyze 
the Boltzmann transport equation. 

In this study, the method of Knudsen iteration, as proposed by 

Enoch (Ref. 1 ) for the hyperthermal regime, is applied to the problem 

of hyperthermal Couette flow with oxidation occuring at the stationary 

surface. The approach utilized involves expressing the free molecular 

distribution function as a delta function, n 6 (v - V ), for the hyper- 

1 — — 00 

velocity particles, and a delta function, n^ 6(v) for those particles which 
leave the stationary surface. Using these approximating forms in the 
relevant collision integrals, Knudsen iteration may be performed with 
the full Boltzmann equation. Closed form expressions are then obtained 
for the heat and mass transfer as functions of Mach number, Knudsen 
number, molecular weight ratios and the relevant physicochemical 
surface parameters. 
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In order to clarify the oxidation processes that occur in rarefied 
flows, the simplest one-dimensional problem has been chosen -- that of 
hypervelocity Couette flow with oxidation at the stationary surface. In 
Figure 1, the flow under consideration is schematically shown. At the 
upper plate (x = d), we consider a source for molecular oxygen, with 
density n^; oxygen is injected through the upper plate and is convected 
to the stationary bottom plate ( x = 0) where it reacts with the graphitic 
surface by means of the following two reactions: 

C(s) + 0 2 - CO z (1) 

c(s) + T °2 "* co (2) 

to form the reactant products CC> 2 an d CO. These products diffuse 
through the streaming oxygen molecules to the upper surface, which 
acts as a sink for the reactant products. The intent of the present 
investigation is to find the heat and mass transfer for a given surface 
temperature, T^, plate velocity, V^, and Knuds en number, Kn^. 

To formulate the model we assume that (i) thermal velocity 
of each component is much less than the upper plate velocity, (ii) 
that the distance between the plates, d, is much larger than all relevant 
mean free paths, (iii) hard sphere intermolecular potential for all 
molecular interactions (iv) diffuse reflection for all species at the 
bottom plate. 

With the above assumptions we can proceed to apply a 
generalized form of "Enoch’s" model to the Couette flow. We first 
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consider the form of the distribution functions to be chosen for each gas. 
For the molecular oxygen (with subscript 3) we write: 

(o)» 


(o) (o)' 

3 

(o)' n 


f 3 = n 3 Ml - ZJ + n 3 6{ v) 


(3) 

(o)" 


where n, = — =— , with, n , the source density of oxygen, and n 

3 2 o 3 

is the free molecular density of unreacted oxygen molecules that leave 

the stationary surface. 

For the reactant products CO (subscript 1) and CO^ (subscript 
2) we write: 


(4) 


(5) 

& £* — 

(o)" (o)” 

where n^ and n^ represent the free molecular densities of reactant 
products that leave the oxidizing, stationary surface. 

Using eqs. (4) and (5) we can now perform Knudsen iteration 
on the Boltzmann equation of each species. This can be symbolically 
written as: 


(o) 

(o)" 


T 

= n l 

6 (v) 

(o) 

(o)" 


2 

= n 2 

6 (v) 


flf. / (o) (o)\ 3 / (o) (o)\ 

4- = jjf, , f . ) + £ J :: ( f. , f. 

i 

j- 1 


ii\ i i / - ji; ij \ 1 ' J 


( 6 ) 


Using the forms for ^ and f^ given in eqs. (4) and (5) the integrations 
specified by (6) are tractable. It can be noted upon inspection that: 
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Using the results of Ref. 1 (Appendix A) we can write for the collision 
integrals : 


/ (o) (or 

J 33 ( f 3 ’ f 3 > 


(o)' (o)" 2 

2n 3 n 3 a 33 

v, V 
3 “ 


6 W 3 ' V“ j X 


( 11 ) 


(o)' (o) 


(l + sgn [V o - v 3 ] ) - n 3 n 3 ir a 33 [ 6(v 3 -VJ+ fi(v 3 )] 


/(o) (o)\ (o) (o)' 

J i3 f i ’ f 3 = a i3 n i n 3 6 


(•■ ■ t) 


( 12 ) 


(o) (o)' 


(l + sgn[V m - v.j) - ir a. 3 n. n 3 V o 6 (v.) ; i = 1 , 2 


/(o) (o,\ 2 (O. (o). / 3 \ 

J 3i f 3 ' f i - a i3 ”l n 3 V 3 V -/ 


(13) 


(o) (o)' 


(l + sgnjv^ - v 3 ]) - 7T a. 3 2 n. n 3 6 ( v 3 - Vj i = 1 . 2 


v. • 

where: , sgn x = n . 

i « 


= 1 ; x > 0 
0 ; x < 0 


(13b) 
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and 2 n a_ represents the cross-section for intermolecular collisions. 
One can therefore write: 


df. 

1 

ix dx 


S i<v) = 



i = 1 , 2 


(14) 


V 


3x 



S 3 (v) = 



It is important to point out that for species 1 and 2 the only 
collisions considered in the model are those between species 1 and 2 
and hypervelocity species 3 particles. In the hyperthermal limit 
( ^ 3 ) this class of collisions will be most important in 

determining the flux of species 1 and 2 leaving the stationary surface 3 . 
However, it must be noted that if the molecular weights, m^ and m^ 
are disparate, there can be a significant difference between the 
velocities of the particles leaving the stationary surface, and collisions 
between these classes of particles can become important. 

3. Free Molecular Solution 

(o) (o) (o)” 

To obtain the values of n^ , n^ and n^ , it is necessary to 
obtain the free molecular solution. As noted earlier, we consider both 
CO and CO^ as being formed at the oxidizing surface. We further define 
reaction probabilities, and y ^ as: 


y 


1 


H 

H 


1 out 
3 in 


H 

H 


2 out 

3 in 
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These rate constants are considered to be functions of surface temperature 
alone, and can be related to the experimental data for graphite oxidation 
at low pressures. To relate y^ and y^ to experiment, we require that 
y^ and y^ reproduce the experimentally determined relation for 
(CO/ CO^) w (Ref. 2), as well as the experimentally determined rate 
relation (Ref. 3). Using a curve fit for the data of Ref. 2, one obtains 
the following relation for y^ and y^ : 



log 10 Z = '( 2>47 / t°r) + 2,1 (18) 


and then by considering the reaction as first order, one obtains a second 
relation between y^ and y^ since: 


n C_ 


► 

o 3 
4 m 3 


*1 + y 2_ 


(19) 


and: 


m = k p 
ox w 3 


( 20 ) 


one obtains: 


3 ( y \ + y 2 


( 21 ) 


Combining eqs. (17) and (21): 


T k w °3 

" 2 = /m 7?~2 

/ 2 \ Z + 1 


m 
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v i ■ 


4 k c _ 

3 w 3 


rn 


m. 


1/2 


T 7T 


m t 


Z + 1 


( 23 ) 


where is the experimentally determined specific reaction rate: 


k = k exp I /RT 

W O V Wj 


(■ E/RT w) 


(24) 


Using the values of y and y given in (22) and (23) one can 

(o)" 

easily make use of the defining relations for y and y to find n , 
(o) (o) 12 3 

n l ’ n 2 : 

_ (°) 
m. C . n. 

H 3in “ “i " — - 1 ;i = 1 * 2 < 25 > 


H 


n 

o 3 


3 in 


(26) 


Therefore: 


(o) 


"i = r i ,-=■ 


n C m 
o 3 3 


2 C . m. 
i l 


/ \ 1 /2 

\> / m 3 \ ; i = 1 

— y i\zr) 


= 1,2 


(27) 


In addition, we have from the conservation of mass: 

H. 


3 out 


H 


3 in 


l - y, - v 2 


(28) 


Therefore: 


(o)" 

n 3 = 2 n Q ( 1 - - y 2 ) (29) 
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4. Results 


(o) (o) (o)" 

With the values of , n^ and determined, one can easily 


write the solution to eqs. (14) and (15): 


X s (v ) 

= -r 1 - 


IX 


+ F. (v ) ; i = 1 , 2 


(30) 


x S (v] 

f 3^— ’ X) = 

3x 


+ F 3 (v) 


(31) 


To determine F^(v), F^(v ) and F^(v) we apply the following boundary 
conditions: 


x = d, f. = 


f 2 « 0 


m (v - V ) 

3 — * oo 

2kT 


r 2ff kT 


TJZ 


; v < 0 
x 


(32) 


x = 0, f. = 


n. 

i 


1 ( 2ffkT ^ 


TJZ 


- mv /2kT 

e ; V > 0 (33) 


and n^, n^, n^ are to be determined by the specification of the inter - 
facial kinetics of oxidation. 

To determine n^ and n^ one employs the definitions of y ^ and y^: 


m. n . C . 

ill 


= y. H . ; i = 1, 2 

i 3 m 


( 34 ) 
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and for using eq. (28): 


C 3 m 3 n 3 


= H 3in (y i + y 2 ' 1} 


( 35 ) 


where: 


H 


3 in 


m f (o, v ) v d v 
3 3 — x 


(36) 


v X < 0 


Using eqs. (34) and (35) one can write for F^(v] 


3/2 n 3 ( 2- V J 


F_ (v ) = n 


3'-' o \ 2frkT 


,3/2 - m, v / 2kT 


'■h 


+ n 


3 \ 2rrkT 


2kT 


1 - sgn v 


1 - sgn v 


(37) 


Substituting the results of eqs. (31) and (37) into (36): 


H 


3 in 


( m 3 "o C 3 / 4 ) 


1 - 


m_ 


t rV n d 

oo o 


3 

, 1/2 


, 1/2 


y i a 31 


y 2 a 32 


(38) 


Using eq, (38), one can then write the expression for the oxidation 
rate with: 


-1 


(Kn ) = d • [ 2rr a, 2 , „ 

\ co / 33 o 


= d • | 2t r n ; S„ = 

o J ® — 
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One can also compute the energy transfer to the surface from the 
definition: 


Q 


3 


f 3 (v, o)m 3 v 3 2 v 3x d\ 3 


v_ < 0 
3x 


Q. 


m^f v i 


v. 

IX 


d 3 v. 

1 


v. < 0 

IX 


i = 1, 2 


Q s Q + 

TOTAL 1 


Q n + CL 


Performing the computation indicated by eqs. (40 - 42): 


Q. 

i 

C n m. V ^ /2 
3 o 1 ® 



CL 


C 3 n o m 3 V » /2 



00 
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C. ENGINEERING CORRELATIONS 
1. Drag 

The drag on an object entering the earth's atmosphere stems 
from essentially different phenomena in the free molecular and continuum 
regimes. In the free molecular regime, the momentum of the hyper - 
velocity free stream particles is transferred directly by impact of the 
particle with the surface. These particle impacts give rise to a 

2 

momentum transfer to the surface per unit time and unit area of P V , 

which can be expressed in terms of a drag coefficient, i. e. , 

C = ^ % 2. However, in the continuum regime, in the 

D 0 V 2 / 2 

oo oo 

vicinity of the stagnation point of a blunt object, the drag on the body is 

2 

caused by the local continuum pressure which is of order: p^ / 2 , 

giving rise to « 1. Therefore, in developing an approximate theory 
to bridge the free molecular and continuum regimes it is necessary to 
consider the drag as a combination of continuum, pressure type drag and 
free molecular, particle impact. Rott and Whitenbury (Ref. 4) have 
succeeded in expressing as: 

C D ’ C D. + „ - C D.) W 

i x F.M. 1/ 

where C represents the continuum value for C and C the 

i F.M. 

free molecular value. In eq. (1), OL represents the fraction of free 
stream particles that arrive at the surface and in general will be a 
function of the Reynolds and Mach numbers. In the free molecular 
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regime when a = 1, one can note that C - C 

8 D d f.m. 

continuum regime a = 0, and - C . Eq. (1) suggests that one 


D. 
1 

C. - 


plot the sphere drag data with 


D F.M. D i 


instead of the con- 


ventional C D . Both Masson et. al (Ref. 5) and Geiger (Ref. f>) have 


C D - C D. 


plotted the data in this fashion I i. e. , 7 "q 


and find the 


D 


F.M. 


D. 

1 


correlation applicable over a wide range of free stream Mach numbers 
and wall temperature ratios. 

Therefore, to develop approximate analytic expressions for 
C D which would be applicable between the free molecular and the 
continuum regimes and be useful over a large range of free stream 


Mach numbers and wall temperature ratios, a curve fit of the data 
presented in Ref. 6 was undertaken. These data are presented in the 
following form: 


C D - C D. 

L vs. Re 

C D - C D. 

F.M. 1 


(2) 


where Re is the Reynolds number behind the shock. Since data is 
sparse for Re g £ 1 0 it has been necessary to use the near free 
molecular theory of Willis (Ref. 7) in the near free molecular regime 
(Re £ . 1) and then interpolate between this theory and the data shown 
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For reference one can write the near free molecular theory of 
Willis in the following form: 

C n 

a 1 - Re [. 835 + 6. 4 S 

C » l ® J 

F.M. 

Equations 3a, 3b, 3c therefore represents a curve fit for the data shown 
in Geiger (Ref. 6), the near free molecule theory of Willis (Ref. 7) 
evaluated at a free stream Mach number of 12 and an interpolation 
between the two. It should be noted that all near free molecular theories 
depend only on the free stream Reynolds number, therefore when these 
theories are plotted vs. Re g there will appear a weak dependence on free 
stream Mach number, so by taking a value of M^ = 12, the result will be 
a good approximation in the range from 7 ^ M^ ^ 15. 


For 10 * Re * 1: 

s 

C - C 
D D. 


F.M. 


D. 


- 6.875 x 10 ^(log in Re 
1 U s 


- 6.875 x 10 log.^Re + 0.84 
& 10 s 


(3a) 


For 1 ^ Re ^ 3: 
s 

C D - C D. 


S - C D 

F.M. i 


- .319(log 10 Re g ) - . 136(log 10 Re s ) + 0.84 

(3b) 


For 3 ^ Re ^ 10 : 


C - C 
D D. 


C D °D 

F.M. i 


7 e 


771(log Re g - . 478) 
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Using a hard sphere intermolecular interaction and considering M^large, 
one can write a relation between Re and Re : 

CO s 


Re 


00 


0. 57 Re M 
s 00 


(4) 


From the experimental data cited in Ref. 3 the value of 

i 

appropriate to sphere drag in hypersonic, high Reynolds number flows is: 


C D = .918 (5a) 

i 

In chosing a value of C , attention must be given to the dependence 

F. M. 

of C on M , and on the surface interaction between the hyper - 

D oo 

F.M. 

velocity free stream particles and the surface. In Ref. 8, Schamberg 
presents a model for the interaction of a hypervelocity particle and a 

surface is given. Using Ref. 8 one can write for C : 

F.M. 


D 


= 2 


F.M. 


1+ ( S co Z 3 )' 1 • 


(5b) 


where 4> is the beam width of the molecules re -emitted from the 
o 

surface, and ^(^ 0 ) is ratio of the axial momentum carried from 
the surface to the momentum that would be carried away if all molecules 
in the beam were aligned with the beam axis, while f(l^) accounts for the 
body shape and reflection law assumed. For diffuse emission one has: 


4 o = 7T/ Z , = -| 

2 

v - « , Uv) = -j- 

4 

C =21+ 

°F.M. L 9y/$~ 



( 6 ) 


(7) 
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To apply eq. 3 to blunt faced bodies ("cigarettes") we assume 

here that the geometry does not have a significant effect on the shape of 

the transition curve between the free molecular and continuum regimes. 

One can, therefore, use values of C and C appropriate to a 

i F.M. 

blunt faced body; these are: 


C 



.918 


C 


d f.m. 


2. 06 


(8) 


In defining Re g in eq. 1, one now uses d, the lateral dimension of the 

blunt face, instead of R_. , the sphere radius. 

o 

It is expected that the correlation equations developed here 
will give a reasonable account of sphere drag and cylinder drag over 
the entire range of free stream Reynolds numbers for 7, and will 

also, to a good approximation, account for wall temperature effects (Fig. 2). 
2. Heat Transfer 

In developing methods for calculating the heat transfer to an 
ablating graphitic surface over the entire range of Knudsen number from 
the free molecular to the continuum regime, one is not only faced with 
the usual difficulties encountered in the transition regime, but 
additionally, one has the complicating physicochemical processes which 
occur at the graphite surface. 

Analyses of heat transfer to an ablating graphitic surface have 
been performed by Scala (Ref. 9) for the high Reynolds number regime 
and in Section (II- B) of the present work the very low Reynolds number 
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the near free and free molecular regimes have been investigated. These 
analyses indicate for the regimes considered the aerodynamic heat 
transfer in the absence of oxidation will give a good representation for 
the heat conducted into the graphite surface in both the rate and diffusion 
controlled regimes. Therefore, we assume here that over the entire 
range of Reynolds numbers the usual aerodynamic heat transfer will 
provide a good approximation to the energy conducted into the surface. 

To obtain a correlation formula for the aerodynamic heat 
transfer we write: 


Q = 


l C H^J L, 
( C H /*■•.)«. 


R. 


R. 


Q 


B. L. 


(9) 


where the ratio: [ C >J R e / (C yj Re represents the 

l ri S'ri.-K. 

ratio of the Stanton numbers in the low and high Reynolds number 

regimes and Q is the boundary layer heat transfer given by 


Scala (Ref. 9) as: 


Q. 


B. L. 


' (it) (33.3*. 0333 [H e -h w 

\ B / air 7 ft. - sec. 


( 10 ) 


with R^ the nose radius in feet, and with given by the modified 
Newtonian distribution as: 


= P, 


1. 55 M 


, 22 cos 0 + p 


sm 0 


36 


UNCLASSIFIED 


01 ) 



i'i RMASSM: 


• • • • t« 

• • • • 

• • • • • 

• • • • 


in units of atmospheres and 0 the angle on the sphere face: 



( 12 ) 


and: 

h = . 3T (Btu. /#); with T in °R. 

w . w w 

air 

= H = . 3T (1 + .2M 2 ); with T in °R. 

ff Q. CO 00 00 

It should be pointed out that the assumptions implicit in the choice of 

[ C H L.R. 

eq. (9) to calculate the heat transfer is that the ratio 


( C H ) H - R - 

is a function of Re alone. It is clear that for hypervelocity flow, the 
s 

above ratio is almost independent of Mach number in the boundary layer 

4 4 2. 

(Re ^10 ) and viscous layer (10 > Re ^ 10 ) regimes; however, in 

s s 

the transition and near free molecule flow regimes there will be a 
dependence. The assumption here is that the dependence is weak. 

To obtain a correlation, expressions are obtained for the 
Stanton number ratio in the various flow regimes and interpolations are 
employed where results are not available. 

In the free and near free molecular flow regimes, the results 
of Willis (Ref. 7) may be applied. One can write, using the results 
of Ref. 7: 


(c u 7Te ) , = /rI [ 1 - . 53 Re 1 

\ H v s / n. f. m. ¥ si ® J 


( 13 ) 
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This result was obtained by a Knudsen iteration of the integral equation 

which results from a modified Krook model. It is interesting to point 

out that one has a dependence on Re in this regime rather than Re , 

® s 

so that to cast eq. (13) as a function of Re g alone requires assuming a 
specific free stream Mach number; for the purposes of this investiga- 
tion we choose: M = 12. 

00 

In the low Reynolds regime we have the theory of Cheng 
(Ref. 10) for a perfect monatomic, non-dissociated gas and that of 
Goldberg and Scala (Ref. 11) for an equilibrium dissociated gas. By 
interpolating between these theories one approximately obtains the 
effects of non -equilibrium dissociation in the gas phase. Further, by 
interpolating between the results of Willis (Ref. 7) and Cheng (Ref. 10), 
one obtains an approximation to the transition regime of rarefied gas 
dynamics. 

One can then write the correlation equations, employing the 
above theories and interpolations (See Figure 3): 


I) 10 -3 i Re £ 10 
s 


C H ' /r % 


L.R. 




7 T exp (* 922 [ lo ®i 0 Re s - l \) 


2. 14 
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(14a) 


II) 10 s Re * 10 
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L.R. 




I 2 

rr I- 405 Nio r 'J +i.«iog 10 R. s+ .5 7 J 
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III) 10 3 ^ Re g 10 5 

C H 2 

, 7=r , L - B - - 272 " I- 055 ('°s,o R %) - • 671 »8 I 0 R % + 4 - 36 1 

(14c) 


It should be noted that the assumptions implicit in using the 

correlation are that M ^7, and that the wall temperature ratio not be 

00 

too large. It is interesting to point out that the use of eqs. (9) and (10) 
in calculating the heat transfer away from the stagnation point implies 
a cos 0 distribution over the entire range of Re^. This is correct in 
both the free molecular and boundary layer limits and can be expected 
to be a reasonable assumption over a wide range of Reynolds numbers. 

For the ’’cigar ette" shaped bodies, one can employ eqs. (9) 
through (14) provided one sets 0 = 0 in eq. (11) and uses the lateral 
dimension of the blunt face rather than in defining Re^. I n addition 
it should be mentioned that for the ’’cigarette” we assume that no energy 
is transferred to the afterbody and that no mass loss takes place from 
the afterbody. 

3. Mass Transfer 

The calculation of the mass transfer from an oxidizing graphite 
surface over the entire range of Reynolds numbers is, in a sense, a far 
more difficult problem than the estimation of either the drag or heat transfer . 
This is because no theories or experiments had existed previously for mass 
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transfer at low Reynolds number, and additionally, the chemical kinetics of 
graphite oxidation has never been experimentally investigated in the 
hypervelocity, rarefied regime. These uncertainties will render the 
estimates that follow approximate, and do indicate the need for more 
extensive theoretical and experimental studies of mass transfer and 
interfacial chemical kinetics over the entire range of Reynolds numbers. 

In discussing the mass transfer problem, it is necessary to 
clearly define the reaction rate and diffusion controlled regimes, since 
either regime may dominate at a given Reynolds number and Mach 
number. 

If the reaction at the surface is proceeding at a slow enough 
rate, compared to the time necessary for convection of the reactant 
species to the surface, one is said to have a reaction controlled process. 

In this limit, the reactant species concentrations in the gas phase are 
generally unaffected by the interfacial kinetics and one can treat the 
surface reactions as being uncoupled from any diffusion processes in 
the gas phase. However, if the reaction rates at the surface are quite 
fast compared to species flow times to the surface, the oxidation pro- 
cess is said to be in a diffusion controlled regime. In this limit dif- 
fusion in the gas phase can be uncoupled from the interfacial kinetics. 

For the oxidation of graphite, Scala (Ref. 9), in the high 
Reynolds limit, has solved the problem of diffusion controlled oxidation 
and has proposed a simple relation for providing a transition between 
the diffusion and reaction controlled limits. Scala obtains for the 
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diffusion controlled oxidation rate: 

T, 1/2 


B. L. 


TTT 


x 6. 2 x 10 


lb 


( 15 ) 


B 


ft 


sec. 


where is given by eq. 11 and for a sphere in units of ft. is the 

nose radius, while for the blunt faced bodies R should be replaced by 

d the lateral dimension of the blunt face. In section (III-B) of the 

present work expressions have been obtained for the diffusion controlled 

oxidation rate in the very low Reynolds number, free molecular and near 

free molecular regimes. In various calculations of oxidation rates 

downstream of the stagnation point performed in this laboratory it was 

found that the mass transfer rates could be correlated by assuming that, 
* 

Q , the effective heat of ablation of graphite was a constant. Therefore, 

in order to construct a correlation for the diffusion controlled mass 

transfer which could bridge the free molecular and continuum regimes, 

* 

it was decided to utilize the same assumption that the Q found by 

Scala (Ref. 9) for the high Reynolds number regime, could be applied 

over the entire range of free stream Reynolds numbers and Mach 
* 


numbers. 


Since Q 
* 


is defined as: 

Q 


( 16 ) 


m D 


one can therefore write: 

Q 


m D = 


❖ 

Q 


(17) 


where Q is the aerodynamic heat transfer rate which may be calculated 
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from eqs. 9 and 14, and Q can be written from Ref. (9) as: 

Q* = 5370 + 5. 37 (He - h . ) BtU /lb 

w air 


(18) 


It should be noted that eqs. (17) and (18) can be put in a form which is 
more amenable to engineering calculations: 

'( C H^ S ) 


m D = m B.L. 


L.R. 




(19) 


where m is given in eq. (15) and the correlation for the Stanton 

B • L# 

number ratio is given eq. (14). 

It is interesting to note that for the free molecular regime 
eq. (19) reduces to the correct free molecular limit, this can be seen 


by noting that: 


S .M. ✓*= 

( C H ^ ), 


, and m. 


2 . 2 


B. L. 


P V 

GO GO 


Re 


s J 


2.2 


so that (m^F.M. 


H.R. 

p V , which is the correct limit. In addition for the 


near free molecular limit as predicted by (19) yields: 


(m) r = p V [l - .53 Re ] 

D n. f. m. oo oo l oo J 


( 20 ) 


which is to within a numerical factor the same result obtained in section 
II-B for the near free molecular oxidation of graphite. Since eq. 19 
correctly reproduces the boundary layer, near free and free molecular 
regimes, it is not unreasonable to assume that it gives a good approxi- 
mation to the diffusion controlled oxidation mass transfer over the entire 
range of Reynolds numbers. 
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In obtaining expressions for the reaction controlled regime one 

encounters large uncertainties which are attributable to the dearth of 

definitive experimental data on the oxidation of graphite when subjected 

to hypervelocity, rarefied flows. 

The existing experimental information on graphite oxidation 

was obtained under gas kinetic conditions (i. e. oxygen molecules reach 

/8kT 

the surface with a thermal velocity C = ^ hrmj ) , these data 
indicate a rate equation of the following form: 


m = k (p n ) 
R w 


( 21 ) 


where — - n ^ 1 and k is a rate constant which is a function of surface 
2 w 

temperature alone. Experimental and theoretical considerations seen 
to indicate that, under gas kinetic conditions: n = , this being 

attributable to a pore mechanism for the oxidation process. During the 
entry of an object into the earth's atmosphere, gas kinetic conditions, 
are generally achieved at the stagnation point after a fully developed 
shock wave has formed. Therefore, for the high Reynolds number 
regime it seems reasonable to utilize eq. 21, with n = ~ > in the rate 
controlled regime, where p^_ 
mentally measured function of temperature. However, as the Reynolds 
number decreases and the oxygen molecules incident at the surface 
retain a significant fraction of their free stream velocities, there may 
be a transition from a pore type mechanism to a surface type mechanism 
which would be first order rather than half order. To clarify these 


- .21 P (eq. 1 1 ) and k is an experi- 
e w 
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questions accurate molecular beam experiments are necessary in the 
hypervelocity rarefied regime. In the light of these uncertainties, we 
choose a rather simple rate equation for the entire range of Reynolds 
numbers, with the understanding that for Re g < 100, the equation will 
become progressively more inaccurate. We write: 

“r ■ V - 2IP / /2 ,22 > 


To provide the transition between the reaction and diffusion con- 
trolled regimes, an addition rule is derived for the boundary layer regime 
and is used over the entire range of Reynolds numbers. We first write 
for the mass transfer: 


B. L. 



(23) 


where the concept of a boundary layer resistance a 

-D # i-j, 

is introduced, and additionally a rate equation: 

“ ■ k »( p ° 2 )° 

' w / 


to mass transfer 


(24) 


Since the mass transfer rate as expressed in (23) and (24) must be equal, 
we equate (2 3) and (24) and obtain the following constraint: 


(i - o m = nc 11 


where: C 


p o, / p o, andTJ 

^w/ 


n-m 

k w a B. L. P 0 2 


(25) 
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Since 0 < C < 1 we may rewrite £ 


as: 


c = 


; /3 * o 


i + a 

Therefore (24) may be rewritten as: 
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m a 




, l7n 

K p, 


+ — 
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Nv p 0 . 


(27) 


l/n 

since k 
w 


• . . 1/m 1/m 

m R and P 0, = m D a B.L. W6Can 


rewrite (27) as: 


(m R ) 


1 , , 1 , 

T7n T : l/n,. TJ 


. ' , m T7r 

(m D ) a B.L. 


(28) 


For the problem of graphite oxidation: m = n = 1/2 and: 

1 


c = 


1 + 77 


(K ^ ) 


Therefore: 


(29) 


m = 


1 1 
+ 


TT72” 


L (m R ) 2 (m D ) 2 J 


(30) 


So that eq. 19 can be used to calculate m jy ec l* 22 can be used in the 

calculation of m and finally eq. 30 determines the total oxidation 
R 


mass transfer from the surface. 
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III. EXPERIMENT 

A. INTRODUCTION 

That phase of the investigation with which the experimentation 
was associated takes no cognizance of the circumstances through 
which particles re-enter the atmosphere. The premise is adopted 
that such particles are re-entering and the investigation seeks to 
examine their behavior under controlled re-entry environment 
conditions. The general solution offered is a machine- computed 
evaluation of the ablation of particles of various sizes following 
selected re-entry trajectories and based on a particular model of the 
ablation mechanism- -that of the oxidation of materials as proposed 
by Scala (Ref. 9)- While the experimental program does not relate 
directly to the computer analysis for re-entering particles, it is 
designed to provide experimental confirmation of the analytical 
approach in two ways. 

First, experiments were performed on ATJ graphite, a well 
characterized graphitic material, to supplement and extend the 
existing data of Diaconis (Ref. 12), which had been obtained in 
support of Scala' s oxidation theory. Data were acquired for the 
reaction rate controlled regime of oxidation for both continuum flow 
and for non- continuum flow. Subsequently, similar mass rate data 
were acquired for a Material "X", furnished by SNPO, and comparisons 
of the performances of the two and with the theory were made by 
empirical evaluation of the test data. 

Second, a study to determine the integrated ablation and the 
shape variation of spherical bodies of different sizes in both a fixed 
attitude and during rotation about an axis perpendicular to the flow was 
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conducted. As a secondary output of the computer program, pre- 
dictions were obtained for the ablation of such spheres subjected to 
environments corresponding to that produced in the arc facility. 

The predicted values were then compared to those actually obtained 
experimentally for ATJ. Tests of several fixed and rotating short 
cylinders constructed from Material n X", were also performed. 

Again, with corresponding spherical specimens, direct comparison 
was made of the results of ATJ Graphite and Material n X M in terms 
of dimensional and mass changes. Further assessment of the com- 
puter analysis was derived through the machine programming of a 
one-dimensional axially oriented "X" cylinder given an input laboratory 
(arc facility) environment. A comparison of the output with corres- 
ponding experimental results was performed. 

B. EXPERIMENTAL EQUIPMENT 
1. Facility Description 

All tests scheduled for the study of ATJ Graphite and Material 
"X" were conducted in a GE-SSL Hypersonic Arc Tunnel. This system 
is comprised of a Tandem Gerdien arc heater, a conical nozzle 
through which the heated test gas is expanded into a continuously 
evacuated test section where the specimens are normally mounted, 
a diffuser aft of the test section, and a high speed mechanical vacuum 
pump. The tunnel is illustrated in Figures 4 and 5. 

The principal element of the test facility combines two Gerdien 
arc units functioning together with a common arc column. The 
arrangement (Figure 6) is such that air is admitted tangentially to 
each Gerdien unit through vortex chambers located between the 
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electrode housings and the plenum chamber. Valving in the air 
system provides a method of dividing the incoming air flows into 
two portions which are subject to close control. The greater 
portion of the air is bled from the system through exhaust ports 
in the electrode housings. In the course of its passage from the 
inlets to exhaust ports, the exhaust air carries with it the carbon 
products produced from the graphite electrodes which would normally 
contaminate the test gas. The remaining portions of the air admitted 
(which constitute the test gas) pass into the plenum from both sides, 
acquiring energy from the arc column enroute. Contaminants, meas- 
ured spectroscopically, constitute less than 100 parts per million 
of the test gas. The rotation induced in the air as it is admitted 
tangentially, is oppositely directed in the two halves of the arc, 
the opposed rotations tending to offset one another as the two flows 
enter the plenum and mix. The combined flow then passes from 
the plenum through the throat of the conical nozzle. 

The nozzle is constructed in two parts, one of which is essentially 
an extension of the other. The up-stream section has an exit diameter 
of 1. 2 n . It is at the exit of this short nozzle that tests at the higher 
model stagnation pressures are performed. The extension, when 
added to the short nozzle, provides a nozzle of 5 n exit diameter. The 
juncture of the two sections is sealed and the internal surfaces form 
a smooth interface in order that the flow will not be disturbed. 

All parts of the arc hardware and the nozzle are cooled with pressur- 
ized high velocity water in order to provide them with maximum 
service life. 
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The test section in which specimens are normally mounted has 
access ports on both sides, top and bottom. Windows mounted in the 
side ports are generally employed for visual observation, optical 
and radiation instrumentation, and in some instances, for motion 
picture camera coverage. The top port is used for camera coverage 
of conventionally mounted specimens. Both top and bottom openings 
are adaptable for auxiliary equipment as needed. 

Normally, specimens are mounted on a water cooled sting which 
extends up-stream into the test section. The sting, operated pneu- 
matically, is capable of holding a specimen in a retracted position, 
out of the flow of heated gases. The controls can be regulated to 
automatically inject sting and specimen into the center of the flow, 
and remove it after a programmed interval. The sting is designed 
to carry instrumentation leads from the specimen through its interior 
to a terminal board outside the tunnel. 

A portion of this program required a special sting to support and 
drive the rotating specimens. This sting was inserted from the top 
port and connected to a motor drive through a gear system (Figure 7). 
Speeds of approximately one and two radians per second were provided 
by interchanging gears. 

Downstream of the test section the system is fitted with a diffuser, 
through which the high velocity of the test flow decreases, with a 
corresponding recovery of pressure. A heat exchanger and a large 
mechanical vacuum pump complete the system. The pump is a 5000 
cfm rotary mechanical booster vacuum pump with two stages, the 
large capacity of which is required to maintain low pressure 
(< 100 4 Hg) in the facility under full flow conditions. 
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2. Facility Operation and Control 

The Tandem Gerdien arc heater is operated on rectified direct 
current from a 500 KW power supply. The power drawn from the 
supply is divided between a ballast resistor bank and the arc column, 
which are connected in series. Power delivered across the arc can 
be varied by changing the resistance of the ballast included in the 
electrical circuit. The resistance of the ballast can be changed 
(between runs) in minimum increments of 0. 065 ohms , to a maximum 
of 1. 300 ohms. Thus, the properties of the test gas can be altered 
in discrete steps as a function of the power reaching the arc. In 
addition, the plenum pressure, and consequently the gas enthalpy, 
can be adjusted by varying the mass flow of air admitted to the arc 
unit and the relative amounts which are permitted to exhaust or to 
enter the plenum. While not customarily done on test programs, it 
has been shown that continuous variation of test gas properties over 
limited ranges is possible by proper control of the mass flows. 

Finally, a discrete change in enthalpy of the test gas can be 
obtained by altering the geometry of the arc plenum chamber or of 
the vortex chambers used in the system. In this program, the 
enthalpy of the test gas was approximately halved by introducing 
vortex chambers with constricting orifices of decreased diameter. 

3. Test Instrumentation and Specimens 

3. 1 Pyrometer: The primary datum acquired on the specimen 

during actual test for this program is a continuous record of surface 
temperature. This is obtained with a recording two-color pyrometer. 
The pyrometer employs two phototubes with maximum responses at 
different wavelengths. One tube has maximum response in the red 
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portion of the spectrum, the other in the blue. The evaluation of the 
relative responses of the two tubes is interpreted as a color temper- 
ature. For models tested in the flow from the 5" diameter nozzle, 
the two color pyrometer views the specimens through a test section 
window. When the 1. 2" diameter nozzle is employed, use is made 
of a special small port built into the arc heater mounting flange. 

3. 2 Total Calorimeter: To measure the enthalpy of the test 

gas, a total calorimeter was used. This device (Figure 8) consists 
of a cylindrical housing containing a labyrinth passage for the test 
gases. The adjacent walls are cooled with pressurized water to 
remove the heat transferred from the test gas. The calorimeter is 
affixed directly to the throat of the nozzle, located in the wall of the 
arc plenum. The heated test gas from the plenum is directed through 
the calorimeter where it transfers the bulk of its heat to the circu- 
lating water. The heat remaining in the gas as it leaves the calorimeter 
is monitored by a thermocouple. Similarly, the temperature rise 
and mass flow of the circulating water are measured continuously 
during a test. Several thermocouples are mounted in the copper 
structure of the calorimeter to evaluate the heat "stored* 1 in the walls. 
The test is continued until all reporting thermocouples equilibrate. 

The evaluation of the reported data, essentially a heat balance between 
the calorimeter elements and the test gas, results in a measure of the 
total heat content of the gas. Since it is estimated that there are some 
minor heat losses which cannot be evaluated, the value of the enthalpy 
obtained will be slightly conservative. 

It should further be noted that the enthalpy determined is a "mean" 
of any distribution of energies at an arbitrary cross section of the 
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nozzle. For example, the heat content of the core gases is not 
distinguished from those of the boundary layer. The latter, being 
adjacent to the cooled nozzle walls, will obviously be lower. In 
addition, other methods, such as spectroscopic evaluation of gas 
enthalpy, indicate that the energy in the core itself is not uniformly 
distributed. Spectroscopic measurements in the past have recorded 
enthalpies about six percent higher than those obtained by total 
calorimetry. 

3. 3 Pressure Probe: As part of the test program it was 

essential to obtain a measure of the model stagnation pressure at 
the various test locations. The principal probe employed to obtain 
these data was a copper blunt cone configuration with a 94 mil central 
orifice (Figure 9). Since under the conditions of the experiment, 
the stagnation pressure is independent of model shape, the same probe 
was used for all measurements. (Other slender conical probes, having 
flat frontal configurations corresponding to the diameters of the mass 
rate specimens and with smaller orifices., were run at the lowest 
density condition to ascertain whether these factors might be import- 
ant to the measurement. Within the required experimental accuracy, 
they were not). The probe was not water cooled. Reliance was 
placed upon the heat sink characteristics of the massive copper body 
to insure its survival in the test environments. In exposure times 
up to 30 seconds, the probe remained intact. 

From the probe, tubing transmitted the pressure to a universal 
manometer rack. The system incorporated an arrangement of 
vacuum valving which permitted both legs of the manometer employed 
to first be evacuated with a vacuum pump. Thus the manometer was 
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referenced to an absolute pressure of less than 10 microns. When 
the pressure probe was inserted into the flow, one leg of the mano- 
meter was valved into the connecting tubing and the pressure was 
read as an essentially absolute value from the displacement of 
the manometer fluid. 

For the lower pressure test condition (5. 0" exit nozzle), the 
manometer working fluid was a vacuum oil of 1. 04 specific gravity 
and very low vapor pressure. For the higher pressure test condition 
mercury was the working fluid. 

3. 4 Specimen Calorimeters : The measurement of heat transfer 

to the specimen configuration, while not essential to the evaluation 
of mass rate data, is considered important in the understanding of 
the behavior of the specimens tested for shape change. This was 
especially true when correlations with the computer analyses were 
sought. 

Calorimetric models corresponding to each of the test specimen 
configurations were exposed to each of their respective test environ- 
ments. Typical calorimeters are shown in Figures 10 and 11. The 
calorimeter bodies were constructed from phenolic nylon, chosen 
for its low thermal conductivity. Into each of these, a copper 
cylinder (in the case of several of the spheres, more than one) was 
inserted. A thermocouple attached to the back face of the copper 
slug recorded the rise in temperature in the copper during exposure. 
By confining the slug within a poor thermal conductor, the heat 
transfer could be evaluated on the basis of a one-dimensional heat 
flow analysis. The one inch cylindrical and spherical specimens 
were equipped with 1/4 M diameter slugs. The 1/4” bodies employed 
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1/8 U diameter slugs. 

To obtain heat transfer to the 1/16" diameter specimen, a rigid 
60 mil copper rod was employed. In order to minimize lateral heat 
input, the rod was coated with a 2 to 3 mil layer of RTV rubber which 
is highly resistant to thermal conduction. The 1 / 1 6 1 1 calorimeter 
appeared to perform extremely well, but in order to provide a 
higher level of confidence in this technique, several calorimeters 
constructed from 1/4" diameter copper rod coated with a correspond- 
ingly thin RTV skin were tested. Results were comparable to the data 
obtained from the l/4 n phenolic nylon cylinders with 1/8 M diameter 
copper inserts. 

3. 5 Test Specimens : Several sizes of cylindrical specimens 

were involved in the acquisition of mass rate data. Sizes included 
diameters of 1/16", 1/4" and 1". For the study of shape change 
characteristics, spherical specimens of one inch and 1/4" diameter 
were included, together with short cylinders having a length to 
diameter ratio (£/d)=2. The configurations and dimensions of the 
specimens are shown in Figure 12. Similar specimens were prepared 
from ATJ Graphite and Material "X" except that only the Hd-Z 
cylinders of "X" were tested. The ATJ Graphite specimens were 
furnished by GE-SSL. A sample of the original block of ATJ material 
is available for property evaluation if required. In general, Material 
"X" was supplied in specimen form by SNPO. 

During the course of the program, various circumstances arose 
which warranted the introduction of several modifications to the 
specimens. Variations were introduced principally to the 1/4" 
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diameter cylinders. It was found that during consecutive tests the 
surface temperatures of similarly constructed models equilibrated 
at temperatures that differed but little from one another. The 
results, therefore, obtained under these circumstances, were 
confined to a relatively narrow range of temperatures. In order to 
extend the temperature range, variations of test specimens retaining 
the basic 1/4" diameter cylindrical shape were constructed (Figure 
13). Objectives sought by making such changes were (a) the alter- 
ation of the thermal conduction paths and (b) a change in the heat 
storage volume of the specimen adjacent to the exposed surface. 

These changes were introduced to the ATJ specimens only. The 
Material "X" specimens were tested in essentially their original 
configurations. While the modified specimens did provide an increase 
in the range of surface temperatures obtained, the effects were not 
as pronounced as expected. However, for the same reasons the 
l/l6 n specimens were also tested in one other configuration (Figure 
13) in addition to the original. In this case, the increase in temp- 
erature range was also relatively small. 

When testing of the rotating spheres was initiated, it was 
immediately evident that the ceramic strut employed did not have 
adequate strength to provide support in the test flow. Consequently, 
a threaded graphite rod was substituted for the ceramic rod 
formerly used. The substitution required that the holes in the 
rotating specimens be deepened and threaded. Reworking of the 
rotating models was necessary for both ATJ and Material "X". 

Both cylindrical specimens and spheres were mounted to the 
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sting assemblies with graphite holders of cylinder- cone configuration. 
Several model and holder assemblies are shown (Figures 14-18). The 
graphite holder was employed to give additional protection to the 
permanent test facility sting. 

3. 6 Motion Pictures: Color motion pictures (16 mm) were 

taken of many of the runs performed at the 5. 0" station. In the place 
of a normal framing speed, time lapse sequences were obtained at 
the rate of one frame per second for some of the longer runs. 
Pictorially, the films may be used to provide qualitative information, 
to verify transient phenomena and to document the tests. However, 
their use for quantitative data reduction purposes is severely limited 
because of the excessive temperature gradients; halation effects 
create the illusion of model enlargement in intensely bright areas. 

A typical effect is the tapered appearance of the two parallel edges 
of cylindrical specimens. In cases where the camera was stopped 
down more, the heated end is more accurately depicted but the 
cooler regions of the specimen do not appear on the film. 

3. 7 General : The facility is instrumented to provide a continuous 

record of arc current and voltage and plenum pressure. • Such inform- 
ation, together with the response curves of the two-color pyrometer 
and calorimeter thermocouple data are recorded on a Midwestern 
multi-channel oscillograph recorder. Mass flows of air are meas- 
ured and monitored by Fischer and Porter float-type flowmeters; 
test flow pitot probe pressure are measured with a manometer system; 
and test cabin pressures are monitored with tilting McCleod gages. 
Pre-set electric timers program operating sequences. These and 
other incidental but pertinent data are hand recorded. 
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C. TEST PROCEDURES 

1. Specimen Location 

From the outset of the program, it was evident that two positions 
would be involved in testing specimens at the lower pressure (5 M 
station) condition. The mass rate specimens, on the coaxial sting 
were adjustable, permitting them to be located at a standard 9/l6 n 
from the nozzle exit. However, the rotating sting assembly required 
that the axis of rotation be located 2. 350" from the nozzle exit. 

Hence the leading faces of rotating one inch spheres were 1. 850 n 
from the exit and those of the 1/4" diameter rotating spheres, 2. 225 M 
from the nozzle exit. 

When testing of the 1/16" diameter mass rate cylinders was begun, 
it became apparent that the for e- shortened elliptical image of the 
stagnation region resulting from the extremely acute viewing angle 
(circa 5°) failed to cover completely the aperture of the two-color 
pyrometer. Since calibration was based on an image which was 
larger than the aperture in the pyrometer, another position (1. 500 M ) 
was established to accomodate the small specimens). 

In the higher pressure (1. 2" station) position all models were 
tested at a location 9/16" from the nozzle exit. The limited field 
of view through the viewing port made this location for the specimens 
mandatory. The field of view permitted a maximum of less than 
1/2" of model erosion to be followed by the two-color pyrometer. 

2. Pyrometry 

One difficulty experienced during testing was a persistant problem 
in aligning the two-color pyrometer to pick up the image of the test 
surface during the initial heating of the specimen. Frequently, it was 
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possible to record a trace only after a manual adjustment of the 
focussing lens, with a consequent loss of as much as sixty seconds 
from the time of model immersion. 

The situation improved when it was discovered that the evacuation 
of the test section produced an appreciable shift in the position of 
the image. The effect was produced by the change in the medium 
through which the radiation passed, with a consequent alteration 
in the refractive angle between the test section and the glass of the 
facility window. When the shift was noted, corrective compensation 
resulted in an improved record of surface temperature from those 
acquired during the early portion of the tests. 

The two-color pyrometer (Ref. 13) used to monitor surface 
temperatures of the test specimens was calibrated against a tungsten 
filament standard lamp. The lamp itself was calibrated with true 
temperature as a function of input current at constant voltage by the 
National Bureau of Standards. This direct method of establishing 
the filament temperature was employed in preference to an optical 
pyrometer calibration in order to avoid the possible human errors 
introduced in reading and in correcting the brightness temperature 
of the filament to a corresponding color temperature. All surface 
temperature data are recorded as color temperatures. In the case 
of ATJ Graphite, the color temperature may be expected to depart 
from true temperature by very little, since the existing emissivity 
data indicate that the spectral distribution approximates that of a 
gray body. While the radiative properties of Material "X" are 
unknown, its composition suggests that surface temperatures 
recorded for it are also good approximations of the true temperature. 

58 


UNCLASSIFIED 


A calibration exercise was necessitated by each major alteration 
of the geometry of the test arrangement. In all, twelve calibrations 
were performed in the course of the program. The responses of the 
two phototubes in the pyrometer were recorded on a Midwestern 
oscilligraph. The deflection of each "pen 11 was measured and the 
ratio of red to blue deflection was calculated. The ratios were 
plotted against their respective temperatures. The resulting graph 
was used to evaluate ratios obtained from the test specimen traces. 

A typical calibration curve is shown in Figure 19. If the source 
radiation increases sufficiently, the tubes will reach saturation and 
the resulting traces will be meaningless. To avoid tube saturation, 
a series of neutral density filters are mounted turret-fashion in 
front of the pyrometer window. As the tube output approaches 
saturation, a filter is indexed to bring the signal down to an 
acceptable level. The calibration technique is, of course, employed 
with each of the filters. 

3. Specimen Heat Transfer 

Measurements of the heat transferred to a cold wall calorimeter 
of the transient type were obtained for specimen configurations 
corresponding to each of those employed for the acquisition of mass 
rate and shape change data. The appropriate calorimetric specimens 
were tested at both the lower and higher pressure test locations and 
in both high and low enthalpy environments. In addition to conven- 
tional stagnation point heat transfer rate determinations on flat 
faced and spherical surfaces, the distribution over a sphere was 
measured. With the 5" diameter exit nozzle, calorimetric measure- 
ments were made at several stations downstream of the exit, since 
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the positions of the several specimens were necessarily different. 

The procurement of a heat transfer rate consists of immersing 
the calorimeter in the flow of hot gas for a very few seconds. Heat 
delivered to the copper insert of the calorimeter by convection is 
conducted the length of the insert to a chr omel-alumel thermocouple, 
the response of which is recorded on the Midwestern Oscilligraph 
Recorder as a function of time. The linear portion of the time- 
temperature trace is evaluated by the method of Sutton and Miksch 
(Ref. 14). Each thermocouple is individually calibrated. 

4. Model Stagnation Pressure 

Stagnation pressure values were required for each specimen 
location and operating condition. The probe, described previously, 
was exposed for a time sufficient to obtain a steady deflection of the 
manometer fluid. In general, less than ten seconds were required for 
the mercury, although oscillations in the vacuum oil were slower 
to dissipate. Test times of up to 30 seconds were used with the oil 
filled manometer. Pressures were measured at each test station 
and for each enthalpy level (Table 1). In the case of the 5 U diameter 
nozzle and high enthalpy operation, values were obtained at several 
distances from the nozzle exit. These are plotted (Figure 20) for 
extrapolation purposes since the probe could not be retracted to 
positions as remote as were several of the spheres. Because simi- 
larity of distribution could be expected and the variation in pressure 
was small, the complete survey was not repeated for the lower 
enthalpy condition. 
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D. EVALUATION OF ENVIRONMENT AND MEASUREMENTS 
1 • Evaluation of Flow Environment 

Data taken from mass rate specimens of several sizes and at 
three operating conditions augmented and extended oxidation rate data 
previously acquired for ATJ graphite and provided comparable infor- 
mation on Material "X". To examine the effect of various flow regimes, 
these specimens were tested in facility environments which, when con- 
sidered with the model sizes, simulated continuum flow including that 
of low Reynolds number as well as near free molecular flow. Differences 
in the rate of mass loss in the various regimes predicted from theoretical 
considerations were sought. 

In order to determine the type of simulation performance which 
the test facility could deliver, it was necessary to define the conditions 
of flow provided at different conditions of operation. Based on facility 
evaluation tests previously conducted in the GE-SSL Hypersonic Arc 
Tunnel (Ref. 15), a series of computations were made to determine the 
Reynolds number at the several test conditions. Although the available 
spectrogr aphic evidence indicates that the flow is in equilibrium as it 
leaves the plenum, sufficient evidence also exists that the flow may 
freeze very rapidly in the vicinity of the throat of the nozzle (Ref. 16). 
Calculation of both the equilibrium and frozen flow expansions through 
the nozzle according to the treatment in Ref. 15 showed that while the 
local free stream properties were quite different at the test locations, 

Llie Reynolds numbers for both frozen and equilibrium flow in the model 
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stagnation region were of approximately the same magnitude, from less 
than 15% for the lower values, to less than 8% for the higher, over the 
pertinent range 500 to 4100 per foot. Consequently, when it was neces- 
sary to introduce a Reynolds number, an average value was used for 
each set of conditions. 

In understanding the character of the test flow environment, it 
is necessary that one be able to recognize the flow regime established. 

An estimate may be made as to whether the gases behave as a continuum 
or as an aggregate of independent particles based on the Knudsen number 
coupled with the Mach number. A knowledge of the Reynolds number is 
extremely useful in evaluating the Knudsen number applicable to the 
various sized specimens. 

Some time ago, Tsien (Ref. 17) proposed an expression for mean 
free path in terms of kinematic viscosity, the specific heat ratio, density 
and sound velocity. X = 1. 255 7) yj y / pu. He further showed, in effect, 
that the ratio of Knudsen number to Mach number was proportional to 
the reciprocal of Reynolds number. By introducing values from kinetic 
theory his expression can rearranged into a convenient form. Kn/M = 
1 . 25 / 7 " /Re. From this equation one may calculate a Knudsen number 
appropriate to either frozen (Kn = 6.5) or equilibrium (Kn = 2.9) flows, 
for the 1/16" diameter specimen. Both numbers are representative of 
conditions generally classified in vacuum practice as a region of transition 
to molecular flow. More recent investigators have advocated the ratio 
of Kn/M as a proper representation of the various flow regimes. The 
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corresponding values are, for Kn = 6. 5, Kn/M = 0.513, and for Kn = 2.9, 
Kn/M=0.604. Hence, it may be concluded that nominally, the same flow 
regime (near free molecular) is predicted regardless of the equilibrium, 
or lack of it, in the heated gas. 

2. Evaluation of Enthalpy 

The nominal magnitudes of total enthalpy of the test gases at the 
two basic operating conditions were measured by total calorimetry. 

Direct measure of enthalpy (h) during each specimen test was impractical, 
however, evaluation of individual test enthalpies was performed in the 
following manner. 

By employing dimensional analysis one obtains the dimensionless 

relationships among the important variables of arc operation. Such an 

analysis is made using electrical power (E^b density (p), mass flow of 

air (m ), nozzle throat area (A), and plenum (total) pressure (p). From 
a 

the analysis the following pi equations are developed: 
hm 


It 


1 


1 T 


2 



”3 = 


Ap 


< E k *.> 


1/2 


Combining 7T ^ and 7T y one obtains 


h = 


TT 
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For this program, where A is constant, the pi factors were eliminated 
by taking ratios with respect to the measured total calorimetric data, 



The subscript m refers to values measured during total calorimetry. 

With the exception of h, all other variables were measured for every 
specimen test, 

3, Evaluation of Data Correlation Factors 

Scala's theory relates mass loss rate and model stagnation pres- 
sure to the characteristic properties of a material as affected by tem- 
perature in the rate controlled regime. Mass loss rate, stagnation 
pressure, and body radius are interrelated in the diffusion controlled 
regime and are relatively independent of temperature. The relation- 

Pb 

ships may be conveniently correlated on a graph of my p vs T^i 

e 

The mass rates per unit area computed by differences must be "normalized" 
with respect to body radius and stagnation pressure in order to be repre- 
sented on such a graph. The normalization of the mass rate data with 
respect to pressure consists of introducing that pressure appropriate 
to a given specimen as determined from the model stagnation pressure 
data. 

The problem of applying an appropriate model radius is a little 
more involved. Geometrically, the flat faced specimen has an infinite 
radius. However, since the body radius enters such a stagnation point 
correlation as representative of the local velocity gradient, one can 
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interpret that parameter in terms of an effective radius. Boison and 
Curtiss (Ref. 18) have made a thorough experimental investigation of 
the stagnation point velocity gradients for model configurations varying 
from hemispheric to flat faced (and cusped). The results of their studies 
relate the velocity gradient (du^/dx) to the geometry of the models. 
Therefore, it is possible to construct from the Boison-Curtiss data a 
graph (Fig# 21) of heat transfer, normalized with respect to that on a 
hemisphere, versus body configuration. Note that for a flat faced 
specimen (x*/R= 0), the heat transfer ratio q^/q ^ ^37. Recalling 

that the heat transfer rate to a body varies inversely with the square 
root of its radius, it is a simple calculation to show the effective radius 
of the flat faced specimen to be 2.46 for R ^ 




Vr 


= JL_ = 0.637 

. 

ff q sph 


assigning q a value of unity. 


R ff " 


1 


(0.637) 


= 2.46 


Simple computations are applicable to flat faced cylinders of different 
sizes. Note that the theoretical analysis in Sect. II C treats with this 
problem by replacing R with the body diameter for flat faced cylinder, 
thus 

R = — R_ 

B f£ (theory) 2.46 (exper. ) 
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Finally, it is necessary to correlate the mass rate data with respect 
to the local flow regime of the model specimens. A portion of the theo- 
retical section of this report deals with the relationships by which the 
correlations are effected. The basic oxidation relationships derived 
by theoretical considerations are constructed for what is nominally 
called high Reynolds number flow. A series of expressions modify the 

fundamental diffusion controlled oxidation theory for various regimes 

-3 +5 

of flow over a range of Reynolds numbers from 10 to 10 . The 

modification is expressed as a Stanton number ratio, designated 

(C TT V^Re ) T „ /(C TT jRe ~ ) _ . For convenience, the functional 

H s L. R. . H s H» R. . 

relationship between Re and (C /rT ) / (C Ae ) has been 

s H s Ijf iv $ H s I !• i\« 

plotted (Fig. 22). The values of the mass loss rates calculated from 
the experimental data reflect the flow conditions at which they were 
obtained. Hence, before they could be compared with theory, it was 
necessary to divide each value by the appropriate Stanton number ratio 
corresponding to the Reynolds number of the test flow in which the data 
were measured. 

The Reynolds numbers per foot determined for the facility opera- 
ting conditions were discussed earlier. Since values of Re^/ft. were 
little different whether computed for frozen flow or for equilibrium flow, 
an intermediate value to both was used in determining the Re g repre- 
sentative of each model size, for each operating condition. 

For oxidation in the rate controlled regime, the data for low 
Reynolds number flow conditions are altered as a function of the effective 
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pressure of oxygen at the model stagnation point. The effective pres- 
sure, treated in Appendix A, can be represented in dimensionless form 
as a ratio to high Reynolds number pressure (Fig. 23) and used as a 
modifying factor to rate controlled mass loss data in a manner analogous 
to the Stanton number factor for diffusion controlled oxidation. In the 
region where transition between reaction rate and diffusion controlled 
oxidation occurred an appropriately weighted mean of the two ratio factors 
was employed. 

4. Evaluation of Mass Loss Data 

The mass loss data acquired in this program was intended to pro- 
vide experimental evaluation of Scala*s theory of graphite oxidation. In 
order to demonstrate the universality of the theory, tests were conducted 
on cylindrical specimens of several sizes under several different con- 
ditions of operation. In all, 89 cylindrical specimens were tested. The 
results of the tests on the various groups of specimens are summarized 
in Tables 2 and 3. 

The useful data obtained from the cylindrical specimens consisted 
of the difference in length resulting from exposure to the test environment. 
The differences were obtained from micrometer measurements taken 
before and after test. Length change at the stagnation point was used 
in preference to measured mass loss because it was impossible to give 
accurate consideration to mass lost from the side wall of the cylinder. 
Typically, the gross mass lost from a cylinder varied between two and 
three times the mass loss on the stagnation region of the specimen, as 
calculated from the length change , using the following expression: 
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Am = pA (Al) 

It is possible that some mass is driven from below the receding 
face of a specimen. Hence a comparative evaluation of mass rate based 
on length change and on mass loss over an area with protected side walls 
would be informative. Provision for such a comparison was made on 
the one inch diameter specimens tested during the program. Since the 
diameter was sufficient, a small insert (1/4 M diameter) was press 
fitted into a centerline cavity in the one inch specimens. Accurate 
mass losses from the inserts were obtained using a semi-automatic 
analytical balance. These are included in the data on Table 4 where 
the ratio of mass loss per unit area based on length change to that cal- 
culated from actual mass loss for each specimen is shown. These re- 
sults are discussed in Sect. Ill, E. 5. 

A few specimens showed a measurable decrease in diameter as 
a result of side wall mass losses. This was especially noticed on the 
1/ 16" diameter models. The net effect of such diametral decrease is 
a tendency for the heat transfer rate at the stagnation point to increase 
slightly (q^ OC 1 /^/rT ). Hence, mass rates for the longer run times 
might be expected to be slightly larger. 

5. Evaluation of the Rate of Mass Loss 

In order to present the results of the tests in a form comparable 
to the theoretical representation it is necessary to obtain the mass rate 
of oxidation for incremental temperatures. While it may be possible 
in theory to evaluate mass rate as a continuous function of suface 
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temperature from motion film records of the tests and time resolved 
surface temperature recordings, the relatively small dimensional 
change in the specimen puts a prohibitive penalty on the accuracy 
attainable. The technique used in this program employed a system of 
length and time differences. 

Specimens were exposed to the test environment for various 
intervals of time, resulting in incremental differences in length change. 
The calculation of rate of mass loss per unit area proceeded from 
taking differences between length losses and times among the various 
specimens inthe following manner: 



The calculated mass loss rate is considered to represent that existing 
at the average temperature for the two test points evaluated, i. e. : 

(T x + T 2 )/2. 

6. Evaluation of Integrated Mass Loss Models 

Specimens of both ATJ graphite and of Material "X M were tested 
to study shape change behavior. Spheres 1/4" and 1" diameter of both 
materials were exposed in a fixed position and at two rotational speeds; 
short cylinders of Material n X" were also examined. The latter were 
studied in two fixed positions, coaxial to and perpendicular to the flow; 
and at two rotational speeds about an axis perpendicular to the flow. 

Evaluation was possible by each of two approaches. First, com- 
parison was possible on a gross weight loss basis. The weight loss from 
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each specimen was determined and recorded (Tables 5-7), and since 
tests were conducted for different lengths of time, an approximation 
of mass loss with time was possible. Second, by studying the changes 
in shape, e.g. (Figs. 24 and 25), one may obtain insight into the cor- 
responding behaviorof re-entering particles. To obtain accurate shape 
change profiles, the ablated specimens were installed on the stage of 
a Jones and Lamson optical comparator. The projected image, magni- 
fied 31.25x was measured to determine the contour to which each shape 
had ablated. On the magnified scale available, dimensions could be 
estimated to values accurate to less than .001" (true scale). 

7. Evaluation of Surface Temperatures 

Many of the surface temperature traces show a "steady-state” 
temperature which departs markedly from a conveniently uniform 
plateau. Variations may be explained by any of three mechanisms. 

It is probable that all three contribute to variations in the trace, and 
presently there appears no way to determine with certainty the primary 
cause of a given fluctuation in the trace. 

First, one may expect the existance of real variations in the intensity 
and spectral quality of energy being radiated from the specimen. Local 
"hot spots" on models have been observed on numerous test runs. These 
intense pinpoints may well cause modification of the spectral distribution 
of the radiation during their relatively brief duration. Such "hot spots" 
are clearly visible in several of the motion pictures. 
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A second source of recorded variation arises fromthe namual 
switching of filters on the pyrometer. Failure to act quickly when a 
build-up in excessive phototube current results in tube saturation and 
produces an interval during which the response curves are not following 
the calibrated values. In some cases, these regions could be deter- 
mined from the original traces, but experience has shown that the 
tubes do indeed recover rapidly and the extra data reduction effort is 
not justified. 

Finally, there is a possible source of error due to the resolution 
of the recording instrument. This source of error is potentially most 
pronounced at the low end of the signal outputs. Thus, when the thick- 
ness of the trace masks an apparently small deflection, an error in 
measurement can affect the reported temperature. 

When the plateau temperature reported for the specimen is highly 
irregular, the best mean value of the temperature was estimated for 
the interval being investigated. In spite of the several sources of error, 
it is estimated that the temperatures used depart no more than 100°K 
from their actual values. 

The validity of measured surface temperature was demonstrated 
dramatically for. one ATJ graphite sphere. A tungsten- rhenium thermo- 
couple was imbedded at the center of a one inch sphere tested at the 
lower enthalpy condition. The thermocouple temperature recorded 
during the entire test was compared with the two-color pyrometer 
temperature measured on the surface. The data plot (Fig. 26) shows 
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a comparatively constant gradient of about 50°K between the surface 
and center of the sphere after the specimen reached a steady state. 
Typical surface temperature traces are found on Figs. 27-29. 

8, Evaluation of Calorimetric Measurements 

Total enthalpy of the test gas was measured for each of the two 
enthalpy levels at which the facility was operated. The results of the 
tests, using the total calorimeter have been summarized (Table 8), 

The evaluation process for data acquired from the calorimeter is a 
heat balance between a heat exchanger (the calorimeter) and the test 
fluid (air) after the system has reached steady state operation. 

The measured value of model heat transfer (Tables 9 and 10) 
are significant because they reflect the changes in model size and 
configuration (Fig. 30), In comparing the measured heat transfer rates 
to an existing empirical correlation equation by Scala (Ref. 20), for 
flight stagnation heat transfer rates, the values measured in the Arc 
Tunnel were found to be significantly lower thanthe values predicted by 
the empirical equation. The difference in heat transfer rates may be 
attributed to differences between the environment of the Arc Tunnel 
and the flight environment for which the correlation equation was derived. 
The equation is based on analyses of the regioncf high Reynolds number 
continuum flow. The flow provided in the tunnel was at low Reynolds 
number. Indeed, facility conditions are such that the boundary layer 
may be fully merged with the shock for all test condition runs, as in- 
ferred by a calculation of theoretical boundary layer thickness and shock 
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detachment distance* This indicated departure from flow conditions 
represented by the heat transfer correlation and exemplified by the 
thick model boundary layer which could produce a reduced enthalpy 
gradient at the wall, probably contributes to the lower (than theory) 
measured heating rates* 

The diminution of heat transfer may be augmented by the frozen 
flow which the tunnel has given evidence of producing* Although the 
model bow wave recompression will affect the gas state, it is not un- 
reasonable to expect a non -equilibrium state in the model shock layer* 
Under these conditions the chemical reactions within the boundary layer 
among the diffusing species could result in incomplete recombination 
at the model calorimeter surface. Further, although the calorimeter 
material (copper) has high catalytic efficiency, minute oxidation of 
the surface could appreciably decrease this efficiency, again resulting 
in a decreased heat of recombination. The measured rates of heat transfer 
were below those predicted by approximately 35%. 

With the exception of the 1/4 U diameter spherical calorimeters, 
all values were plotted directly as calculated from the recorded traces 
by the method of Sutton and Miksch (Ref* 14). Aside from experimental 
errors, it has been ascertained that the values for the spherical one inch 
specimens are conservative by 7%. In attempting to determine apparent 
discrepancies in the heat transfer to the 1/4" spherical calorimeters, 
attention was directed to the size of the calorimeter insert with respect 
to the dimensions of the sphere. A 1/8" diameter slug was installed in 
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a 1/4" sphere, with its face contoured to form the continuous spherical 
surface. Since it represents such a large portion of the entire exposed 
hemisphere it must also sense local heat transfer rates which are less 
than those at the stagnation point. Hence the response of the thermo- 
couple represents an average heat rate less than that incident at the 
stagnation point. The distribution of heat transfer rates on a spherical 
surface is well known. Indeed, as part of this program, the distribution 
on a one inch sphere was measured (Fig. 31). The mean heat rate re- 
corded by the calorimeter was employed in the determination of the 
maximum, or true stagnation heat transfer rate. It was found that the 
experimental distribution of the heat transfer could be approximated 
reasonably by a Gaussian type expression: 


q o e 


-Pitt 


which, when doubly integrated gave the average heat transfer rate: 
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Using information from an experimental study of heat transfer to a 
hemispherical surface (Ref. Zl) the constant P (dimension = L *) was 
evaluated. The resulting equation was then used to calculate a stagnation 
point heat transfer rate from the "average value” reported by the calori- 
meter. The analysis showed that the average value was 53. 5% that of 
stagnation. Similarly, on the one inch sphere, the average value was 
93% that of stagnation with a 1/4” diameter slug insert. 
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The profile of heat transfer to a sphere (Fig. 31) was obtained 

by using several spherical calorimeter bodies tested on separate runs. 

o 

The first has 1/4" diameter inserts installed at 0 (stagnation point), 

60° and 120°. Others included 30°, 45° and 90°. The results of the tests, 
when combined, formed the complete profile. Note that several positions 
were tested twice and repeated the local value with fidelity. Considered 
collectively, the heat transfer measurements show some scattering 
which may be attributed to calibration error, thermocouple attachment 
variations, and trace readout error. 

In spite of a reasonable agreement of heat rates obtained with the 
RTV silicone rubber coated 1/16 M copper calorimeters with the empirical 
curve, serious doubts were raised as to the effectiveness of the RTV in 
inhibiting lateral heat input. As a check of this, several 1/4" copper 
inserts were coated with the silicone rubber and tested. Their re- 
ported heat transfer rates were comparable to those measured with 
calorimeters of 1/4" diameter phenolic nylon containing 1/8" diameter 
copper inserts, thus confirming the validity of the technique. 

As a consequence of ablation, the spherical oxidation specimens 
gradually altered in shape. The departure from a sphere was most 
pronounced on the fixed specimens. To determine what effect the 
shape change had produced on the incident heat transfer rate, calorimetric 
bodies were made in the shape of the ablated spheres. A stagnation point 
calorimeter was installed and measurements indicated that on bodies 
which had suffered the maximum erosion, the effective heat transfer 
rate had decreased to the equivalent of a flat faced specimen. 
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The calorimeter bodies duplicating the ablated specimen were 
prepared by molding. First a plaster shell was formed around the 
ablated graphite sphere. In turn, a casting of RTV silicone rubber 
was made using the plaster shell. Before the rubber was poured, a 
dummy insert with connecting rod was affixed to the mold to provide 
for the later insertion of the copper insert and thermocouple. The 
casts faithfully reproduced all surface irregularities, so that the 
RTV calorimeter body was an exact replica of the original ablated 
graphite specimen. 

E. DISCUSSION OF RESULTS 
1. Oxidation of ATJ Graphite 

Once the numerous mass rate values had been computed (Tables 
11 - 16), they were plotted as a function of reciprocal temperature 
with the objective of determining the activation energy (E) and the re- 
action rate coefficient (k^) for the equation. 

• V /G" -E/RT 
m = k JP e 
o v e 

For pure reaction rate controlled oxidation, the given equation can be 
rewritten in the form of the integrated Arrhenius equation, which re- 
quires that the natural log of m ' vs 1/T shall be a straight line. * 
Examination of the data plotted in this form (Fig. 32) revealed that they 
did not follow the predicted form of the Arrhenius equation. 

1. In this report, direct reference to unadjusted mass rate data 
employs the symbol m' ; corrected data is denoted by the con- 
ventional m . 
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Since it was fully expected that some of the specimens would have 
been oxidized under conditions of diffusion control, a departure from a 
linear plot was not surprising. The deviation obtained tended to confirm 
the presence of the diffusion controlled activity. Coupled with the scatter 
of the data points, the narrow range of mass rates for which linearity 
might be assumed made it difficult to plot a straight line which could be 
considered properly representative of the reaction. In lieu of an estimated 
slope of great uncertainty it was decided that a value would be chosen which 
was typical of graphite materials. The slope is a measure of the activation 
energy. Consultation with the available literature led to a choice of 44. 0 
kilocalories as a good representative value. A straight line of this slope 
was constructed to pass through the bulk of data points estimated to be pre- 
dominantly reaction rate controlled. Based on the constructed line, the 
reaction coefficient was determined: 

k Q = 4. 42 x 105 lb. / sec. ft. ^ atm. ^ ^ 

The equation for the reaction rate oxidation of graphite thus determined, 
was combined with the theoretical diffusion controlled oxidation expression 
using the transition equation derived in the theoretical portion (Sec. II C. 
eq. 30). Curves were computed for each type specimen and test condition 
from the combined expressions (see Appendix B) and were overplotted on 
the data points (Figure 32). It was immediately evident that in the diffusion 
controlled regime, the data points were high with respect to the theory. 
With the knowledge that the test facility provided environments which 
were not characterized by high Reynolds number flows, these differences 
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were expected. 

The theory developed by Scala (Ref. 9) to describe the oxidation of 
graphite was based on flow of high Reynolds number . In this report, 
that theory is extended to regions of low Reynolds number and to mole- 
cular flow. The extension provides that data of both rate controlled 
and diffusion controlled oxidation must be adjusted in order to compare 
it with high Reynolds number continuum based theory, or vise versa. 

In their final forms, the modifications consist of applying a Stanton 
number ratio correction factor to data in the diffusion controlled regime 
or an effective pressure ratio correction factor to data in the reaction 
rate controlled regime. Both these factors are functionally dependent 
upon the Reynolds numbers corresponding to the test flow conditions. 
From calculated nominal values of Reynolds number per foot, corres- 
ponding to the various test conditions, the shock Reynolds number for 
each model was determined. The appropriate correction factors for 
either oxidation process were then obtained from the graph of Stanton 
number ratio vs. Re g (Figure 22) or the graph of | p^ / .21 P ot J 
vs. Re g (Figure 23). Data point values must be divided by the cor- 
rection factors determined from these graphs. 

In examining the mass rate plotted against 1/T, it may be seen that 
many of the points lie near the transition region of the theory curves. 
For these points, correction factors are needed which are intermediate 
between those appropriate for diffusion control and for rate control. 

To accomodate these points, the theoretical expressions were examined 
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and the contribution of each factor was weighted for each point in the 
transition region. Refer to Appendix C. Once the appropriate correc- 
tion factors were determined for all regions containing data, they were 
applied, each to its corresponding datum point, to determine the mass 
rates in terms of Scala's original theory (Tables 11-16). 

A second graph of m vs. 1/T was prepared showing the adjusted 
data overplotted on the theoretical curves (Figure 33 and Appendix B, 
eq. 4). It may be seen that the diffusion controlled data are decreased 
and conform reasonably well to the predicted performances. In general, 
the rate data are increased in value, although the changes in position 
with respect to the theory curves are less pronounced because of the 
near -vertical slope of the reaction rate portions of the curves. 

It is interesting to note that when the data are corrected to correspond 
to the theory, the shift of reaction rate data may be sufficient to warrant 
the construction of a new, displaced, slope for the activation energy and 
alteration of the weighting factors in the transition region. Extrapolated 
to its extremity, this becomes an iterative process which ceases when 
the required degree of accuracy has been attained. In the present case, 
the scatter of the data points is sufficient to make further approximations 
unjustified. 

Having shown that the mass loss rate data does indicate reasonable 
concurrence with the theory, the final step was to prepare a plot of the 
data normalized with respect to pressure and specimen radius. How- 
ever, it was evident that data corrected to correspond to the theory 
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curves would tend toward a more crowded pattern. In order to better 
examine the relative correspondence of theory and experiment for the 
individual test conditions, the theory curves were modified as the cor- 
rection factors dictated, while the original mass rate data were norma- 
lized with resepct to pressure and body radius only, (Tables 11-16). 

The results were plotted ( m' = J Rg/ P e vs. T w ; Figure 34) to pro- 
vide a clearer picture of the degree to which theory and experiment 
agree (Appendix B, eq. 7). 

Finally, the appropriate correction factor was applied to each datum 

point normalized with respect to pressure and specimen radius (Tables 

11-16), and the computed results were incorporated into a plot 

( m - J R /p vs. T ; Figure 35) in which the theory curves remain 
' v B e w 

in the form calculated from the unmodified continuum theory (Appendix 
B, eq. 6). The theoretical curves predict a dependency of data on the 
body size only in the reaction rate controlled regime but independence 
with respect to both pressure and body size in the diffusion controlled 
regime. The methods employed in reducing the data are described in 
detail in Appendix C. 

Aside from the scatter of data it may be seen that there is indeed a 
trend which places the 1/4" specimen data between the majority of one 
inch and of 1/16” specimen data. In addition, the data for all three model 
sizes are in approximately the proper relation to one another as predicted 

o -j / 2 1/2 

by the theory. Note also that the 6. 2 x 10 lb. / ft. sec. atm 

plateau is placed in the midst of the data points corresponding to diffusion 
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The experimental results assuredly concur with the theory when a 
mean value of the data is considered. The difficulty of controlling pre- 
cisely the operating conditions during the course of a run produces 
appreciable scatter. Hence the exact reproduction of the theory curve 
by experimental evidence is exceedingly difficult, although continued 
experience with this type of experiment has led to improvements in 
performance. 

2. Integrated Mass Loss - ATJ Graphite 

To obtain integrated mass loss data and shape change information, 
tests were made on spherical specimens, (a) fixed with respect to the 
flow and (b) rotating at two different speeds about an axis normal to the 
flow. Data were obtained from measurement of the profiles of specimens 
exposed for incrementally varied time intervals and from determinations 
of the gross loss of mass experienced by each. Profiles were obtained 
from measurements of the images projected by an optical comparator. 

Profiles of the several initially one inch diameter fixed spheres were 
constructed to scale (Figure 36). From these , a plot of dimensional 
change vs. central angle was constructed for each of the increments of 
test time (Figure 37). The central angles were measured with the stag- 
nation point at the zero angle position. It was also convenient to con- 
struct a plot of dimensional change vs. time for particular central angles 
(Figure 38). Time was used as a parameter for this graph. 

Examination of the length loss vs. central angle plot indicates that 
the removal of mass from the specimen reflects the variation of heat 
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transfer to the spherical surface (Figure 31). It may be recalled that 
the theoretical portion of this study postulates a constant effective heat 
of ablation (q*). The experimental data obtained on the spherical 
specimens suggest that this is a valid assumption, since, if q* is to be 
constant, the rate of heat transfer experienced at a point must be a con- 
stant ratio with respect to the mass lost at that point. Of course, as 
the model changes shape the level of heat flux will also change, de- 
creasing as the spherical surface becomes flattened and increasing as 
it becomes reduced in size. 

Profiles were constructed for the 1 M diameter spheres rotated at 2 

rad. 

radians per second ( q) = 2 sec. Figure 39), for which sufficient 
time variation in the tests provided adequate data to define trends in 
performance. Three intervals of test time were available. The time 
dependent graph of length change (Figure 40) does not exhibit the double 
inflection that was observed on the stationary spheres, within the 
accuracy of the measurements. One may infer that the effective heat 
transfer rate on the stagnation region, being considerably reduced 
because model rotation distributes the accepted heat more equitably, 
prevents the increase in length loss from progressing fast enough to 
clearly define a change in slope. A tendency of the slope to decrease 
at the longer test times reflects the gradual flattening of the equatorial 
region of the sphere. The cross plot of length loss vs. central angle 
(Figure 41) shows the equatorial flattening more clearly. The final 
shape assumed by the rotating specimens was that of a prolate spheroid. 
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The progression in shape with central angle again relates to the heat 
transfer curve, as was noted for the fixed spheres. 

Similar data were acquired for the 1/4 M diameter spheres, both 
stationary and rotating. Complete profiles of the stationary specimens 
(Figure 42) were constructed,' and the plots of length change vs. time 
(Figure 43) and of length change vs. central angle (Figure 44) exhibit 
characteristics similar to those of the one inch spheres. One 1/4" 
diameter stationary specimen (not included on the preceding graphs) 
was tested for 398 seconds, resulting in virtual destruction in that time. 
Only a slender annular piece remained, attached to the ceramic model 
support. 

When profile measurements of the 1/4” diameter spheres at two 
rotational speeds were made, it was found that the dimensional changes 
were essentially identical. Hence, the profile measurements were 
averaged and a composite profile was constructed (Figure 45). Despite 
the similarity in observed dimensional change, it is interesting to note 
that positive differences in mass loss were detectable. These will be 
discussed later. With only two increments of test time available, the 
length of the 1/4" rotating spheres were not plotted as a function of time. 
However, a plot of the measure of length vs. central angle (Figure 46) 
shows that length loss decreases as one progresses from the equator 

toward the pole ( Q = 0° >0 = 90°), similar to that which was 

observed on the one inch spheres, but to a lesser degree. 

Finally, a graph of integrated (gross) mass loss from the various 
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spherical bodies was prepared (Figure 47). Note that for clarity, the 
mass loss scale used for the one quarter inch specimen data is 16 times 
that used for the one inch specimen data in that their values were in- 
creased in inverse proportion to the l/4 n diameter spherical surface 
area compared to that of the one inch spheres. In several instances 
there were insufficient data points to define the shape of a curve. In 
these cases, the illustrated curve was constructed in a shape analogous 
to those which were well defined. 

It may be seen that for both size models, there is a definite and 
consistent trend for the mass loss to decrease as the speed of rotation 
increases, although at a diminished rate at the higher speeds of spin. 

This suggests a method of approximating the degree to which spin rate 
will affect the effective rate of heat transfer. It is evident that the 
effect may be somewhat variable with time, since the bodies ablate to 
different shapes when fixed than when spinning. Assuming the plotted 
mass losses to be directly proportional to the heat transfer at any time 
(q* = constant), the apparent reduction in q (spin factor) was deter- 
mined (Figure 48) for the two spin rates employed. The value of q on 
the fixed sphere was taken as unity. 

Note that the ratios obtained tend to increase with time. This is 
obviously due to the decreased heat transfer rate (the reference q 
of our ratio) evolved on the fixed sphere, as it flattens during ablation. 
Later in time, the more slowly ablating rotators will also develop large 
flattened equatorial zones and hence the slope may be expected to decrease. 
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Given time for the fixed sphere to be ablated to beyond the original center, 
the slope would probably reverse, since the transverse diameter of the 
fixed specimen would start to decrease. 

If one assumes the ratios may be extrapolated to zero time (dotted 
lines), the relative value df the spinning sphere to the fixed sphere, 
before the latter has suffered appreciable shape change, may be 
estimated. These values were obtained for both size spheres, and 
a cross plot of the ratios vs. the model rotation speeds (Figure 49) 
shows that the ratios quickly approach constant values at little more 
than co = 4 radians / second. 

There is yet another method of considering the effect of rotation. 

That is, one may look at the local mass loss per unit of area as re- 
flected by local dimensional changes. At any instant after time zero, 
the stagnation point of the fixed sphere will have receded more than 
the stagnation "point" (of these, there are an infinite number) on the 
rotating sphere. To examine this method of presentation, the stagna- 
tion point length loss data of the fixed and the rotating spheres were 
plotted (Figure 50). Ratios of the length lost by the rotating sphere 
to the length lost by the fixed sphere for each size specimen were 
computed and plotted as a function of time (Figure 51). There was 
insufficient evidence to separate effects due to differing rotational 
speeds but the time dependence of the ratios were again evident. The 
similarity in shape of the two curves is fortuitous; considering the 
scatter of the experimental data, differently shaped curves would have 
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been much more probable. Nevertheless, one may attribute the time 
dependency of the ratios to the different deformation experienced by 
rotating and fixed spheres. 

It appears that sphere size is also a factor. The same conclusion 
may be drawn from the ratios obtained for gross losses in mass 
(Figure 48). Note that both methods, of presenting the results, mass 
and stagnation length loss, indicate a higher ratio for the smaller 
spheres for long exposure times, and both indicate a "cross-over" in 
the ratios for 1/4" and 1" spheres at about 250-300 seconds. 

Either method can be used to help predict performance of a rotating 
sphere of a given size, if one has information on the performance of a 
stationary one. The former permits estimation of the net loss in mass 
and the latter, estimation of the equatorial dimensional change. 
Obviously, either should be deducible from the other, but the connective 
relationships are not easily evaluable. 

3. Oxidation of Material "X" 

After the individual test results had been paired and the mass loss 
rates determined (Tables 17-21), an Arrhenius plot of m vsl/T was 
constructed (Figure 52) which exhibited a pattern whereby the data 
fell within a narrow, essentially linear band. There was no positive 
evidence that diffusion controlled oxidation had developed. 

Assuming the data to be primarily in the reaction rate controlled 
regime, a slope was constructed through the data points and the 
reaction rate coefficient (ko) and the activation energy (E) were 
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determined. The reaction equation derived from the experimental 
data was 

n / e — 1 9700 /T 

m - 2090 V P e 

An approximate analysis (Appendix D) provides conservative 
prediction of the plateau established when diffusion processes control 
the oxidation: 


V R -3 3/2 1/2 

B ^8.5 x 10 lb/ ft sec-atm 

P e 

This and the empirical rate controlled equation were combined using 
the transition equation (Section II C, eq'n 30) to provide a means of 
calculating the theoretical oxidation curves for Material "X" (See 
Appendix B). These were overplotted on the experimental data 
points (Figure 52) where it may be seen that in general, the data did 
not depart markedly (beyond experimental scatter) from the "theory" 
curves. 

However, since the data were acquired under conditions other than 
high Reynolds numbers, the values obtained required adjustment in 
order to be represented as equivalent high Reynolds number mass 
losses for proper comparison to the theory. The required modifications 
were made (Tables 17-21) and with the exception of the 1/4" models 
in the short nozzle, all data were corrected as predominantly rate 
controlled. Hence the appropriate effective pressure ratio corrections 
(Figure 23), as dictated by the corresponding shock Reynolds numbers, 


87 


UNCLASSIFIED 




MHUI A QQinri* 

--'"'Ointa 


• • 



were applied. Since the 1/4 M specimens in the short nozzle appeared 
to fall in a transition region, the weighting factors applied depended 
upon the extent to which each process, rate control and diffusion 
control, contributed during the transition. The factors were applied 
to the mass rate data and the resulting modified data were plotted 
on a second graph of mass rate vs. reciprocal surface temperature, 
(Figure 53). It may be seen that the corrections applied did not 
produce a major shift in the rate data with respect to the theory curves. 
Again it may be emphasized that when a higher degree of accuracy is 
justifiable, the reaction rate curve may be iterated and the individual 
data points adjusted by new correction factors consistent with the 
altered theoretical curves. As with the ATJ graphite, no iteration was 
applied to the Material "X M data. 

The mass rate data were normalized with respect to pressure and 
model size and the corrective values were applied to the normalized 
data (Tables 17-21). The final comparison was a graph of the norm- 
alized mass rate data plotted against the corresponding temperatures 
(Figure 54). The empirical-theory curves are shown on the same 
graph. It may be seen that the data correspond reasonably well with 
the "theory” curves through the reaction rate regime. There is some 
disagreement between the two with regard to transition to the diffusion 

controlled regime. A number of data points appear to project the 

-3 3/2 1/2 

reaction rate to values greater than the 8.5 x 10 lb/ft sec-atm 
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conservatively predicted for diffusion controlled act ivity (Appendix D). 
Indeed, the experimental data does not appear to indicate the incipience 
of diffusion controlled activity up to the limit of recorded values. 

4, Integrated Mass Loss - Material M X M 

Spherical specimens of Material "X" analogous to those of ATJ 
graphite were tested under corresponding conditions. In addition, 
several short cylindrical specimens were exposed to the simulated 
re-entry environment. Fixed cylinders were oriented both coaxial 
to and perpendicular to the flow. As with the spheres, the tumbling 
cylinders were tested at two rates of rotation (end over end). The 
axis of rotation was perpendicular to the flow. Figure 55 shows the 
full variety of "X" specimens, one of each type (including the mass 
rate specimens) employed during the experimental program. 


From profiles obtained using the optical comparator (Figure 56), 
the length losses of the one inch fixed spheres at two enthalpies were 
determined and plotted as a function of central angle (Figure 57). Both 
sets of curves exhibit the double inflections, indicating that the length 
changes are dependent upon the distribution of heat transfer across 
the surface of a specimen. With only two increments of time, an 
accurate representation of the time dependent variation of length loss is 
not assured; however, for illustrative purposes, curves are drawn 
(Figure 58) showing trends analogous to those determined for the ATJ 
graphite with respect to the stagnation point length changes. There is 
no reason to believe that the general shape of the time dependent 
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behavior of Material "X" differs drastically from that of the ATJ. 

Since so little change in radial length was observed for the one 
inch rotating spheres of Material "X", even after nearly 800 seconds 
of exposure, the enlarged profiles illustrated on Figure 59 show only 
the maximum change as compared to the original contour. The actual 
incremental radial length losses for the rotating spheres, based on 
profile measurements are plotted in Figure 60. A comparison with 
the stationary sphere (run at the corresponding flow conditions) is 
made on Figure 6l for the stagnation point location. The large 
difference in radial length loss which shows up in this comparison does 
not also apply to the integrated mass loss since a great deal more 
surface is affected on the rotating spheres (Figure 67). 

The 1/4" diameter spheres, both fixed and rotating, were also 
measured and profiles were prepared (Figures 62 and 63). Based on the 
profiles, graphs were constructed depicting radial length loss as a 
function of central angle (Figures 64 and 65). When the time dependency 


of stagnation point length change for fixed and rotating 1/4” spheres 
was plotted (Figure 66), it was necessary to construct approximations 
of the curve shapes through the limited number of data points. Again, 
it is evident that stagnation point length change comparisons do not 
represent those for integrated mass loss (Figure 67), although certain 
correlations can be made. 

To summarize, the radial length measurements are subject to some 
error because of the irregularities of the ablated surfaces of Material ”X” 
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(Figure 68). In addition, the time resolved length change curves, 
constructed through two or at most, three data points are not 
necessarily representative -- they are contoured in a form analogous 
to those which were obtained with the ATJ specimens. Thus the 
question of whether the integrated mass loss data obtained on the 
spherical specimens tend to confirm any of the observations with 
regard to the radial length changes is one which must be considered 
within the limitations of the relatively small number of data points 
obtained. Total (integrated) mass loss, plotted as a function of time 
(Figure 67) for both size spheres -- fixed and rotating, depicts the 
actual mass losses of the l/4 n spheres multiplied by the ratio A' /A 
(surface area of a 1" sphere divided by that of a 1/4" sphere) to 
allow presentation of both sets of data to a comparable scale. Note 
that the normalized mass losses are approximately two to one, which 
are proportional to the heat transfer rates applicable to the two 
different sizes of spheres (Table 10). They also correlate well 
when compared on the plot of rn\j R b / P e vs ^w {^S ure 54), the 
steady state surface temperatures averaging about 200°K higher for 
the 1/4" spheres (R^-, 125/12 ft. ) in the region which is clearly 
reaction rate controlled. 

In similar manner to the comparison plot of q^ / vs time (Figure 48), 
constructed for. ATJ graphite spheres where q* = q/m is considered to 
be constant, the ratios of A m r / A m f ^ or Material "X" spheres, taken 
directly, are plotted against time to show that the trends are 
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significantly different and that the q* for "X" is probably not a constant 
for these conditions (Figure 69). 


Finally the results of the short cylinders ( / d-2) were plotted in 

terms of mass loss versus time (Figure 70). In terms of the maximum 
time of test, they may be ranked according to the mass loss which 
each incurred. The cylinders placed in a fixed position with axis 
perpendicular to the flow lost the most mass. The next greatest mass 
loss was sustained by the cylinder tumbling at to = 1 rad/sec. followed 
by that having an to = 2 rad/sec. , entirely consistent with all previous 
comparisons. Least total mass was lost from a cylinder mounted 
coaxial to the flow ( M cigarette 11 ). Again, based on the ATJ spheres 
for illustrative purposes, the cylinders were assumed to have doubly 
inflected slopes. The graph suggests that mass losses from the 
tumbling cylinders may be compared to those of stationary cylinders 
with axis perpendicular to the flow in a manner completely analogous 
to the spheres, as discus sed previously. 

It is interesting to observe that the coaxially mounted cylinder 
lost the least mass. Clearly, in spite of the relatively high heat rate 
one may associate with the one quarter inch diameter exposed face, 
the net mass loss was diminished because only a relatively small 
surface area was exposed to the higher heat rate. 

To assess the mass loss contributed by the side walls of the 

\ 

coaxial cylinder, the net change along the axis was measured and 
converted to an equivalent mass for the volume lost assuming no 
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diameter shrinkage. It is perhaps surprising to note that about 64% 
of the mass lost comes from the side wall. This evidence is sufficient 
to confirm that an idealized one dimensional heat conduction solution 
predicting material performance for the "cigarette" case would be 
entirely inadequate. 

5. Comparison of Material "X" and ATJ Graphite 

The density of ATJ graphite is 108 lb/ft^; that of Material "X" 
was determined to be 134 lb/ft^. The evaluation was based on the 
average of fifteen measurements of weights and volumes of the 1/4" 
diameter inserts supplied for the test program. These values were 
used in the reduction of all data. 

An obvious inquiry might be whether the rates of mass lost from 

test specimens of the two materials was proportional to their densities. 

Inspection of the graph (Figure 71), prepared to illustrate the normalized 

mass rates of both materials, indicates that such was not the case. The 

lower density ATJ lost more mass per unit time than did Material "X" 

below about 1800°K. The empirical rate equation derived for ATJ 

e -44000/RT 

graphite, m = 4.42 x 10~ 5 e » can be compared to the rep- 

resentative "slow" and "fast" reaction graphites referred to by Scala. 

(Ref. 9). The reaction rate obtained for ATJ falls between them, 
although it is nearer in value to the "slow" graphite. The activation 
energy of a material determines the characteristics of the exponential 
portion of the rate expression. Note that the activation energies of the 
two materials are similar. Their slopes on an Arrhenius plot would 
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be nearly parallel. Hence the reaction rate coefficients are predominant 
in determining whether they are relatively ’’fast" or "slow” in reacting. 
Clearly, the Material H X M is very slow. Indeed, it has a rate approaching 
that of pyrolytic graphite. For similar activation energies, the rate is 
slower as the reaction rate coefficient decreases. 

The ATJ graphite develops a well defined region of diffusion con- 
trolled oxidation. As temperature increases above approximately 1500°K, 
a value of mass rate is established which is independent of temperature. 
The experimental values obtained for Material M X' r do not give positive 
evidence of leveling off within the range of data obtained. Consequently 
the exponential rise in mass rate of Material "X" crosses the constant, 
diffusion controlled value for ATJ. Thereafter, as the temperature 
increases, mass is lost from Material "X" at a rate greater than from 
ATJ. It is expected that Material "X n might conservatively exhibit a 
diffusion controlled region, as shown by the curves constructed for 
m”^/^Rg/P e = 8. 5 x 1CT 3 lb/ft^^ sec. -atm^^ (Appendix D). Due to 
an unforeseen shortage of the smaller mass rate specimens of "X" 
which developed as a result of the significant difference in oxidation 
performance experienced between "X" and ATJ, the experiments with 
U X" were not extended far enough to establish a defined value for dif- 
fusion control. Perhaps, in the future this should be done, along with 
a thorough analytical evaluation of the chemistry of reactions. 

For both materials, a comparison of mass lost from the 1 /4 ,f 
diameter inserts of one inch mass rate specimens was made with their 
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measured changes in length. To put the length and weight losses on a 
common basis they were both converted to weight lost per unit of area. 
The results were tabulated (Table 4). It may be seen that the ratio of 
unit mass loss-based on length, divided by unit mass loss-based on 
weight, consistently resultedin values of less than one for the ATJ. 

This result suggests that there was a small amount of mass driven 
from below the surface of the graphite. In appearance, the texture 
of the graphite surface became coarsened as a result of exposure to 
the arc heated air tending to suggest development of a porous surface 
of decreased density. 

The same effect, to a much more pronounced degree, was 
observed on the surface of Material ’’X’ 1 . Indeed, on the comparator 
enlargements, the profiles of Material "X" specimens were so ir- 
regular that frequently the mean surface was poorly defined. Because 
of the gross porosity of the inserts employed in the one inch diameter 
Material U X M mass rate specimens, there was a greater scatter in the 
ratios. Nevertheless, the mean value of the ratios falls below that of 
ATJ graphite. The relative effect on Material M X' r is about 1-1/2 times 
that on ATJ graphite. However, when one compares the general texture 
of the two surfaces (Figure 68), the severe irregularity of Material ,f X ,r 
suggests that penetration below the mean surface was effected at more 
localized sites , whereas the penetration of the ATJ graphite appeared 
to be completely uniform in distribution. 
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With respect to the integrated mass loss specimens (1 /4 ,f and 1" 
diameter spheres) of each material, both stationary ( 60 = 0) and spinning 
{OJ = 1&2 rad/sec), comparisons were made of the losses in radial 
length and weight sustained by the ATJ graphite and by Material "X". 
First, the total mass loss equivalent to the volume removed from each 
specimen was estimated by assuming simple geometries to describe the 
volume changes. For the stationary sphere, the difference between the 
volume of a hemisphere and one half the volume of an ellipsoid of revolu- 
tion was used. For the rotating sphere, a uniform spherical shell of one 
half the equatorial radius change in thickness was employed. Using the 
densities for the ATJ graphite and Material "X", the corresponding total 
mass losses were determined. They compared within a few percent of 
those actually measured by weighing the spheres before and after test. 

The radial length changes sustained by the one inch diameter 
stationary "X" spheres were taken from the data on Figure 58 and 
ratioed to corresponding ATJ values from Figure 37. The ratios 
increase from 0.36 to 0.65 as a function of time to 900 seconds. Thus, 
the radial length changes of ATJ graphite are always greater than those 
of Material "X", but the relative loss becomes less with increasing 
time. On performing such an exercise, one can see that these values 
are consistent and check with the mass rate curves (Figure 71). The 
increase in ratio can be attributed to two factors: 

(1) The shape change (flattening) of ATJ spheres is more 
severe. Consequently a decrease in heat transfer to the surface 
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(2) As test time is increased the ATJ graphite performs longer 
in the diffusion controlled regime while "X" continues to perform 
along the rising rate controlled curve. The observed temperatures 
(^1700°K for ATJ; ^ 1820°K for H X n ) are such that the mass rate 
curves do not cross. Therefore, consistent with the actual measure- 
ments, the ratios appear to approach unity. 

A similar examination of the two materials can be conducted for the 
rotating spheres. For example, when considering the one inch diameter 
spheres, the ratio varies from 0.32 to 0.22 as time increases from 200 
to 900 seconds, (Figures 61 and 40). Casual inspection might suggest 
that the results of the spinning spheres are inconsistent with those ob- 
tained for the stationary spheres. However, it should be recalled that 
Material "X" displays a surge of lost material during the very early 
part of each test. Thus, the ratio should be expected to decrease as 
the mass rate from Material "X u falls to a more uniform, u steady - 
state" value. Note that the initial surge of lost material is confirmed 
by the mass loss data vs. time curves for Material "X" (Figure 67). 

Study of the results obtained from the one quarter inch diameter 
spheres, although based on a more limited number of test points, 
indicates that the trends observed with the larger spheres persist. 

On comparing the data for both the stationary and the rotating spheres 
(Figures 43, 46 and 66), one finds that for the longer exposure times, 

M X M looses more radial length than ATJ. This is consistent, since 
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the smaller spheres experienced surface temperatures which are higher 
than the cross over temperature found on the comparison plot of the mass 
rate parameter (Figure 71). 

Total mass losses were also compared. For the stationary spheres, 
an increase of ratio (mass loss of "X"/mass loss of ATJ) was again noted. 
The ratio ranges from 0.36 at ZOO seconds to 0.74 at 900 seconds. Hence, 
the results are approximately the same as those based on the radial length 
change ratios. 

However, in the case of the rotating spheres, the weight loss surge 
was relatively more predominant. Consequently (for the one inch diameter) 
the ratios gave a maximum at approximately 400 seconds and thereafter 
declined. 

While the ATJ graphite and Material M X U did not ablate in identical 
modes, in general each material is self consistent when considered with 
respect to weight loss and dimension change. The peculiarity of Material 
M X M in exhibiting a large loss of material for a short period during the 
early portion of test exposure is interpreted to be the principal cause of 
the differences in performance between the two materials, exclusive of 
those variations which are associated with their oxidation mass rate 
behavior as illustrated in Figure 71. 
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IV. SYSTEMS STUDY 

A. INTRODUCTION 

Appropriate flight trajectory, heat transfer, heat conduction and 
mass transfer equations have been incorporated into an IBM 7090 
computer program which enables the re-entry behavior for specific 
shapes and materials of interest in this study to be predicted. The 
effects of variable drag coefficient, fragment weight and shape change 
on the flight trajectory have been included. Two specific fragment 
shapes may be analyzed one-dimensionally: 

(1) a rotating sphere whose actual behavior is assumed to be 
simulated by a monotonically shrinking sphere experiencing an 
average heat and mass transfer which is some fraction of the 
stagnation value uniformly distributed over the spherical 
surface, and 

(2) a cylinder whose axis is aligned with the local velocity 
vector (zero angle of attack) experiencing a uniform heat and 
mass transfer on its front surface which is a spherical segment 
whose radius of curvature equals the diameter of the cylinder 

(or such other radius as may be more suited to a speaific analysis) 

and remains constant while the length of the cylinder diminished 
during re-entry. 

Data input required include the following: 

(1) re-entry velocity, path angle and altitude. 

(2) thermal conductivity, specific heat as a function of 
temperature and density. 

(3) activation energy and specific reactivity. 

(4) selection of spherical or rectangular coordinates (sphere or 
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cylinder shape). 

(5) initial dimensions of sphere and/or cylinder. 

(6) appropriate averaging coefficient, <I>, which modifies 
stagnation heat and mass transfer. 

Data output includes: 

(1) instantaneous temperature distribution through the material. 

(2) instantaneous mass loss rate, heat transfer rate, and weight 
of the body. 

(3) instantaneous radius of sphere and length of the cylinder. 

(4) atmospheric properties such as stagnation pressure and 
free stream Reynolds number. 

In addition to the one-dimensional idealizations, a multi-dimen- 
sional heat conduction solution involving a moving boundary was 
programmed for IBM 7094 digital computation. This program enabled 
a better approximation to be made regarding the shape changes 
associated with a re-entering sphere and also provided a tool for 
analytically simulating laboratory shape-change experiments. The 
program was written in generalized coordinates to allow flexibility in 
selecting coordinate systems for analyses. The theoretical heat and 
mass transfer relationships developed in Section II are incorporated 
as boundary conditions. 

B. TRAJECTORY CONSIDERATIONS 

A Two degree of freedom trajectory was selected to represent 
the motion of fragments re-entering the atmosphere. The effects of 
lift forces, rotating earth and winds were neglected. Lagrange's 
equations of motion were used to determine the differential equations 
for a point mass in polar coordinates. Figure 72 is a schematic of the 
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motion of a re-entry fragment. Lagrange's equations state that the 
rate of change of the generalized momentum is not in general equal 
to the generalized force, but equal to the force plus the derivative 
of the kinetic energy with respect to the corresponding coordinates. 
Thus, the Lagrange equations 



k = r, © 


mr ^ 0 + 2mrr 0 = F^ 
in polar coordinates where 


T + 1/2 mV 2 = 1/2 m [r 2 + (r 6) 2 ] 


ST 

af 


mr 


ST 
d r 


mr 


ST _ 

sT - 


2 • 

mr 0 


S T 


o 


(1) 

( 2 ) 

< 3) 

(4) 


F is the total force measured along the radius vector 

\2 


F = D sin y - mg 
r o 




F = 
r 


V C n A sin y - mg 


m 


The negative sign fixes the convention as downward negative. 

V = [r 2 + (r e > 2 ] 1/2 

r r 0 

Sin y = - — ; cos y = -y- 

Combining equations (2), (5), (6) and (7) yields 

5,2 1 - r.2, , X,Z- ] l/z( C D A ) . / R e\ 2 

r-re = - 7 P„ Lx +{r 0 ) ] ^ 6() r - j 


(5) 


(6) 

(7) 
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Fg is the moment about the center of the earth. 

F q = - D c° s y ( r ) = - I p 9 V 2 C d A cos y (r) (9) 

The negative sign in equation (9) fixes the convention as positive to the 
right. Substituting into equation (3) yields 

r 0 + 2r 6 = j p 9 [r 2 + (r &) 2 ] 1/2 g Q (10) 

Equations (8) and (10) have been programmed with appropriate initial 
conditions as input to define the trajectories used in the systems study. 
The effects of decreased mass and frontal area as well as variable drag 
coefficient as described in Section IIC. 1 are continuously determined 
throughout the trajectory. 

C. TRANSIENT HEAT CONDUCTION 
1. One Dimensional Idealization 

The conduction of heat throughout an arbitrary volume may be 
written in vector notation as 


pc P If- = V *( KVT ) an 

where K is the thermal conductivity coefficient which may be both 
position and temperature dependent. Figures 7 3 and 74 illustrate the 
monotonically shrinking sphere and one dimensional cylinder ideali- 
zations which were analyzed in the systems study. 

For the monotomically shrinking sphere, equation (11) becomes 
in polar coordinates 

pc P T T = ^2" 7TC (f R2 !r) (12) 

and the following boundary and initial conditions applied: 

T (R, 0) = T e 
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I t), t] 


Q w (s,t) - ae w t w 
av 


(0, t) - 0 (13) 

It should be noted that the symbol s = s(R, t) denotes the instantaneous 


surface position which moves with a velocity 


R = 


W 


av 


(14) 


W 


and the symbol (s,t) represents the suitably averaged heat 

av 

transferred to the sphere surface and rh^ the suitably averaged 

av 

rate of mass loss. 

For the one dimensional cylinder, equation (11) becomes in 
rectilinear coordinates 
ST 


Pc p 'St 




(15) 


Q w (s.t) 

av 


CT€ w T w 


and the following boundary and initial conditions applied: 

T (x, 0) = T q 

ts(x,t),t] = - j^( 

Ci.t] = o 

In like manner, the symbol s = s(x, t) denotes the instantaneous 
position of the surface which moves with a velocity 

m w 

• av 


( 16 ) 


% 

and the symbol (s, t) represents the suitably averaged heat 


(17) 


av 


transferred to the stagnation region of the cylinder. 

Equations (12) and (15) have been included in an IBM 7094 
digital computer program together with trajectory equations developed 
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T (R,0 ) = T 

o 

Cs (R,6),t] =i [Q (0, t) - a e T 4 ] 

3 R K w w w 

lim T = 0 < T < «> , 0 < 6 <77 
R -0 

where Q (0,t) is the local heat transferred to the sphere surface and 
w 

the temperature at the center of the sphere is bounded. 

Appendix E contains comparisons of exact and numerical solutions 
for both equations (12) and (18) utilizing boundary and initial conditions 
(13) and (20). 

D. OXIDATION OF SPHERES 

Table 22 summarizes the results of the systems study for the 

idealized monotonically shrinking sphere. Figures 77 through 91 present 

the results of the parametric study graphically. The material thermal 

properties which were used in the study are illustrated in Figures 92-94. 

It should be noted that two sets of chemical kinetic data were used in the 

2 i/2 

parametric study, viz., k^=2090 and 2.0 lbs/ft -sec-atm , and E=31900 

and 18000 cal/mole . These represent the results of correlations from 

experimental data obtained during the final months of the study and the 

results of an earlier experimental correlation, respectively. The diffusion 

-3 -3 3/2 

controlled regime constants of 8.5 x 10 and 6. 2 x 1 0 lbs /ft -sec-atm 
also represent final and earlier approximations for Material "X" (See 
Section II E). For ATJ graphite, the "slow” kinetics (Ref. 9) were 
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in Section IV. B. and heat, mass loss, and drag coefficient correlations 
developed in Section II. C. 

2. Multi-dimensional Considerations 

A more exact analysis of the mass loss associated with objects 
re-entering the earth* s atmosphere requires a multi-dimensional 
formulation of the heat conduction within the object. A multi-dimen- 
sional heat conduction subroutine was developed in generalized 
coordinates to enable several geometric shapes to be analyzed with 
only minor modifications required in the program. 


Figures 7 5 and 7 6 illustrate the rotational modes of the sphere 
analyzed multi-dimensionally. For an isotropic sphere, equation 
(11) becomes 



boundary and initial conditions apply for the sphere rotating about an axis 
o 

of 90 to the flow direction: 

T (R, 9 ) = T 

o 

H: Cs(R.8 > ).t]=i [Q w (e.t) -a e w T w 4 ] 

av 

§^- ( ftt) = 0 , (19) 

where Q (0 , t) is the suitably averaged heat transferred to the sphere 
av 

surface. For the case where the axis of rotation is aligned to the flow 
direction (or stationary), the following boundary and initial conditions 
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o 


The average heat transfer and mass loss rate was assumed to be 37. 5% 
of the stagnation value as discussed in Appendix F. 

Figure 77 illustrates the percent mass loss for the monotonically 
shrinking sphere as a function of initial ballistic parameter (or size). 

It should be noted that the parameter, W/C^A, is equal to 4 yO w R/3C^ 
for a sphere and equal to for the cylinder, and is a more 

convenient parameter for describing the particle performance character- 
istics. The particle sizes investigated ranged from 1/1 6 to 1.0 inches 
in diameter. Figures 78 and 79 illustrate the instantaneous radii of the 
re-entering spheres for both the earlier and final experimental values 
of reaction-rate and diffusion controlled constants for Material "X", 

For the range of particle sizes investigated the presence of a critical 
boundary of particles having an acceptable initial size was not found; 
only slight mass loss was experienced by the Material "X" spheres. A 
maximum mass loss of approximately 14% was found for a sphere having 
an initial radius of 1 / 4 in. and a re-entry angle of Z. 5 . Since the percent 
mass loss, defined as the ratio of the change in particle weight to the 
initial weight, is inversely proportional to the ballistic parameter and 
must be zero when the ballistic parameter is zero, the location of a maximum 

2. The kinetics determined experimentally for ATJ (Section III E) a short 
time after the systems studies had been started, differ only slightly 
from the "slow” kinetic values assigned by Scala, thus an adjustment 
was deemed unwarrented. 
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percent mass loss is not surprising. A comparison of the instantan- 
eous radii of a 1/4 in. ATJ and a Material "X” sphere re-entering 
is illustrated in Figure 80. Figures 81 and 82 illustrate the effect of 
varying the initial re-entry path angle, and also, the averaging factor 
which modifies the heat transfer and mass loss rate to the monoton- 
ically shrinking sphere utilizing the earlier Material "X" reaction-rate 
and assumed (graphite based) diffusion-controlled regime data. The 
percent mass loss increases with increasing angle of re-entry. Since 
the level of mass loss for this material is low, the fact that the actual 
averaging factor was unknown introduced no major uncertainty. Figures 
83 thru 86 illustrate the trajectory parameters, surface temperatures 
and mass loss rates associated with the re-entry of a typical ATJ 
and Material n X" sphere compared in Figure 80. Figures 87 and 88 
summarize the heat transfer data for the spherical re-entry objects. 

In Figure 87, the presence of a maximum in surface temperature 
illustrates a transition between very small and very large particles 
(corresponding to low and high values of the ballistic parameter) 
from a heat conduction point of view. For extreme values of the ball- 
istic parameter, very small objects have negligible heat capacity 
while the very large particles act as heat sinks. Since the heat transfer 
rate increases with the ballistic parameter, the presence of a maximum 
surface temperature is anticipated. The results of a study conducted in 
Reference 24 found a similar maximum in surface temperature for re- 
entering stainless steel discs. From Figure 88 it may be concluded 
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(1) The maximum heat transfer rate increases with increasing 
re-entry path angle for a given value of the ballistic parameter, 


w/c d a. 

(2) The integrated heat transfer rate decreases as the re-entry 
angle increases. 

(3) The maximum and integrated heat transfer rate increase 
with increasing value of the ballistic parameter. 

Figure 89 illustrates the maximum temperature gradients experienced 

during the re-entry of the spherical and cylindrical objects. (Note: 

Cylinders are discussed in Section IV-E). The small temperature 

gradients experienced by the spheres suggest that a simpler approach, 

viz. , assuming infinite thermal conductivity or equivalently an isothermal 

sphere, would have given reasonable values of mass loss for values 

2 

of the ballistic parameter less than 4 lbs/ft . Figure 90 illustrates 
the effect of varying the ballistic parameter on a velocity-altitude -time 
plot. The effect of varying the initial re-entry path angle for a constant 
initial ballistic parameter is illustrated in Figure 91. 

E. OXIDATION OF CYLINDERS 

Figure 74 illustrates the idealized cylinder analyzed in the 
systems study. Table 23 summarizes the results of a parametric study 
on the behavior of re-entering cylinders. The average heat transfer and 
mass loss rate was assumed to be equal to the stagnation value assoc- 
iated with a radius of curvature equal to the cylinder diameter. (See 
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Appendix F). Figure 95 illustrates the percent mass loss as a function 
of the ballistic parameter. It should be noted that the percent mass 
loss is inversely proportional to the stagnation region radius of 
curvature. A maximum percent mass loss of 17% was found for a 
1/4 in. dia. cylinder of Material "X" having a length equal to 1/2 in. 
using the final experimental values of the reaction-rate and diffusion- 
controlled regime constants. This value is 20% higher than the maximum 
value found for the monotonically shrinking sphere. Figures 96 and 97 
illustrate the instantaneous length of the re-entering cylinders for both 
the earlier and final experimental values of the reaction-rate and 
diffusion-controlled regime constants for Material "X", Figures 98 
and 99 illustrate the trajectory parameters , surface temperature and 
mass loss rate for the cylinder experiencing the maximum percent mass 
loss. Figures 100 and 101 summarize the heat transfer correlations 
for the earlier cylinders investigated. It should be noted that the 
surface temperature levels are considerably higher than those of the 
sphere for approximately the same value of ballistic parameter. For 
larger cylinders, the possibility of transition to the sublimation regime, 
although not investigated in this study, could possibly occur and result 
in considerably higher mass loss. 

F. OXIDATION OF SPHERES (MULTIDIMENSIONAL APPROACH) 

The development of the multidimensional, moving boundary 
computer program during the latter months of the study enabled an 
assessment to be made regarding the conservative or non-conservative 
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nature of the monotonically shrinking sphere idealization used to 
conduct the systems study. Two rotational modes, illustrated in 
Figures 75 and 76 were investigated, viz. , a sphere whose rotational 
axis is aligned with the flow (or equivalently a stationary axisymetric 
sphere) and a sphere whose spin axis is aligned 90 to the direction of 
flow. These two orientations were assumed to represent the maximum 
and minimum mass loss situations associated with a re-entering sphere. 
A comparison of the three re-entry modes for an initial 1/2 in. radius, 
Material "X" sphere follows: 


(A) 


(B) 


(C) 


Mode 


Monotonically Sphere, spin Sphere, 

shrinking sphere axis aligned spin axis 

(random move- with flow 90°to flow 

ment of spin axis) (= stationary) 


Percent Mass Loss 

12.4 

10. 2 

3. 2 

Maximum Surface Temp- 
erature (Stag. pt. ) 

4100° R 

4500° R 

3580° 

Averaging Factor for 
Heat & Mass Transfer 1 

r m w 

V X£ - av =. 375 

w 

o 

o 

« 

II 

s 

m 

= — 

m 

Change in radius at 
stag. pt. 

. 022 in. 

.056 in. 

. 008 i 


av 
■ ( 9 ) 


= . 375 
aligned 


Figure 102 illustrates the surface temperature variation at the stagnation 
point for each of the above cases as a function of re-entry time. The 
average heat and mass transfer factor for case (C) of . 375 was verified 
to some extent for the stagnation point during the analytical simulation 
of laboratory shape -change experiments (See Section V). 
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It may be concluded from the above comparison that the integrated 
mass loss obtained for the monotonically shrinking sphere, also utilizing 
an averaging factor of . 375, is non-conservative since it predicts a 
mass loss even greater than the maximum mass loss associated with 
the multidimensional analysis, (Case B). Thus, the percent mass losses 
obtained in the systems studies conducted in Section IV -D are non- 
conservative and should be adjusted downward by using a lower averaging 
factor in order to produce percent mass losses which are bounded by 
cases (B) and (C). 
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V. COMPARISON: THEORY-ANALYSIS WITH EXPERIMENT 

In the process of attempting to compare computer evaluated pre- 
dictions of material performance in the facility environment with 
actual experimental results, evidence was acquired that theoretical 
expressions describing the flight environment were not directly 
applicable to the conditions of the test environment. In addition, 
it was demonstrated that the heat transfer relationships introduced 
into the computer program were not well suited to the prediction of 
actual shape change on spherical specimens exposed in an arc 
facility. Recognition of the differences between the test and flight 
environments led to a modified comparison of the loss in length of 
Material "X" cylinders as determined by computer calculation and 
by experiment in which the analytical expressions were altered to 
provide a heat transfer rate similar to that measured experimentally. 
Thereupon the results produced were found to compare quite 
favorably. While the findings of the comparison should not detract 
from the validity of the theory as applied to re-entry, it demonstrates 
the need for the use of actual measured stagnation heat transfer rates 
in the computations describing performance of materials in a facility 
test. 

The first computer comparisons were made for stationary (or 
coaxially spinning) ATJ spheres since their behavior was more 
completely documented experimentally. For these calculations. 
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only the experimentally measured gas enthalpy and model stagnation 

pressure were introduced into the computer program. All other 

conditions, such as heat transfer and pressure distributions, were 

evaluated from the programmed theoretical expressions. The 

predicted performance of the stationary spheres as obtained from 

the computer data is illustrated (Figs. 106-108). Performance data 

were obtained for 600 seconds of exposure. On comparison with the 

experimental data (Fig. 109) it was found that predicted radial 

length changes were significantly larger than those measured except 

for the region near the stagnation point. Refer to Figs. 38 and 107. 

In addition, surface temperatures were well above those recorded 

during test (T £2 3100°R at 9 = 0 ). 

w exper. 

A similar disparity existed between computation and experiment 
for the rotating ATJ spheres. These were programmed for calcula- 
tion with the same data as was provided for the stationary spheres. 
One additional factor was introduced. Since the spinning of the 
sphere reduces the effectiveness of the heat transfer to the specimen, 
it was necessary to introduce into the computer program a reduced 
effective heat transfer rate. The choice of the reduced rate was 
dictated by consideration of both analytical evaluations (Ref. 26) and 
by evidence from the experiment. A factor of 0. 375 based on 
stagnation point information was used. The predicted performance 
of the spinning sphere obtained by computer is illustrated 
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(Figures 110-112). The run was evaluated for 600 seconds. As in 

the case of the stationary sphere, it was found that predicted mass 

losses were high compared to those actually measured as inferred 

by AR measurements (Figure 113). Refer also to Figures 41 and 

111. Similarly, surface temperature values exceeded those 

measured (T C& 2900°R at 0 = 0°). It is most significant 

exper. 

however, that the relative mass loss curves for both fixed and 
rotating spheres are similar. This consistency of results suggests 
that the spin factor introduced for the rotating spheres was of 
approximately the correct magnitude to predict satisfactorily the 
effect of rotation. 

When one attempts to analyze the differences between prediction 
and performance, the greater radial length losses obtained on the 
computer immediately recall a similar discrepancy between 
measured and predicted heat transfer rates, discussed in Section 
III D. 8. There, it was found that Scala's correlation equation, 
described as showing good agreement with experimental flight data, 
predicted heat transfer rates significantly higher than those measured 
in the test facility. The distributed heat transfer to a sphere has been 
plotted (Figure 114) showing both the predicted and the measured 
values to illustrate the differences. It may be seen that, relatively, 
the heat transfer rates exhibit characteristics which would produce a 
ratio curve similar to the mass loss ratio curves as defined by radial 



114 


UNCLASSIFIED 





• • 

• • 
• • 

* c 

• • 


:WU$WJEi 

• • • • • • • • 

• * • • * • V • 


• •• 

• • • 

• • • • 

• • • 

• • • 


length loss. Note that if the measured values of heat transfer were 
introduced into the computer program to reduce the mass rates pre- 
dicted, the prediction would be much more representative of the 
results actually obtained. 

The variation in the relative magnitudes of the heat transfer rates 
as a function of angle with stagnation point suggests the only obvious 
area in which the flight theory may not be entirely satisfactory for 
application to the facility experiment. The forward hemispheric 
region of the specimens are, according to analytic representation, 
subject to a heat transfer distribution based on a modified Newtonian 
pressure distribution. Hence the variation in heat transfer rate is 
approximately proportional to the cosine of the angle made with the 
stagnation point. Experimental data suggests that a distribution 
related to the cosine square of the angle might be more appropriate 
for the test facility. This tendency seems to be confirmed by the 
increase through a maximum of the relative values of the heat 
transfer as one progresses around the sphere toward the 90°position 
(Figure 115). In summary, a better computer prediction of facility 
performance can be obtained by incorporating a representation of the 
measured heat transfer rates into a program. 
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experiment, the body-cylinder acted as an oxidation guard ring, also 
providing thermal balance, for the 1/4" diameter cylindrical inserts, 
from which the measured length changes were obtained. Hence, the 
heat transfer through the insert approximated the one dimensional 
situation. Since heat transfer is a variable dependent on the nose 
radius of a specimen in the computer program, an effective radius 
was chosen for the calculation which led to a rate of heat transfer 


approximating that to which the experimental specimens were sub- 
jected. In turn, the measurements taken from the inserts were 
considered to be representative of results obtained from an actual 
specimen, the overall diametral dimension of which was one inch. 

This is, of course, realistic, since the incident flow distinguishes 
the specimen as a whole. The other computer input data, gas enthalpy 
and stagnation pressure, correspond to the conditions of test. 
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A graph (Figure 117) was prepared on which both computed and 
experimental length changes have been plotted. A mean curve was 
passed through the data points and extrapolated linearly to 600 
seconds. It is immediately evident that the computed values show 
better correspondence with the experimental results, although the 
experimental data indicates a change in length early in time which 
the computer prediction does not. 

Hence, as a result of the comparisons, it may be concluded that 
the analytical evaluation of the oxidation process presented here, 
coupled with material properties and trajectory information, may be 
incorporated into a computer program to yield reliable predictions 
of actual re-entry behavior. The differences existing between flight 
environment and the test facility environment require that inputs to 
the computer program be sufficiently modified to compensate for the 
discrepancies between flight and facility environments, particularly 
in the case of heat transfer. 

Of particular interest is the apparent departure of the mass loss 
profile plots relating to rotating spheres of Material "X" from a 
correspondence to the heat transfer distribution such as occured on 
similar rotating spheres of ATJ graphite. Refer to Figures 31,39,41, 
45,46,59,60,63 and 65. Note that for "X" the radial length changes 
are relatively uniform over most of the spherical surface. This is 
not so for ATJ, and may well be attributable to the differences in 
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thermal properties between the materials. By implication then, for 
bodies which do not constantly present the same face to stagnation 
heating, materials which behave the rmo-chemically like Material "X" 
appear to lend themselves to re-entry degradation solution by the 
simpler one dimensional trajectory analysis - -that of the monotonically 
shrinking sphe re - -pro vided, of course, that the proper heat transfer 
averaging factor ( M ) is applied. ATJ graphite, for example, and 
perhaps many materials, require for analyses the more elaborate 
multidimensional solution. Stationary bodies (or equivalently those 
spinning about an axis coaxial to the flight direction) should be 
treated multidimensionally regardless of the constituent material. 
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Throughout the text of this report the various considerations have 
been discussed, and where possible, conclusions drawn in context with 
the discussions. The primary conclusions are, however, summarized 
as follows: 

(1) Over the range of re-entry situations considered, and probably 
for any shallow re-entry trajectory where the velocity is orbital or less, 
small objects of Material M X M will not decrease via oxidation mass 
transfer to an acceptable impact size with respect to the potential 
radiological hazard. Thus for this material, a critical boundary of 
minimum initial particle size is not defined and the computation of 
residual activity inventories based upon integrating over a statis- 
tical distribution of re-entering particles would be currently academic. 

(2) Upon experimental determination and/or verification of 

E and C for a specific graphitic material, a mathematical expression 
describing the aerothermochemical mass transfer characteristics 
of that material can be derived and applied in the convenient form 
m V r b /p e vs. T^. Further, although Reynolds number and surface 
pressure have an appreciable effect on the oxidation behavior in the 
non-continuum flow regimes, the results can be expressed and 
compared in terms of continuum oxidation theory; conversely, the theory 
can be used to define the oxidation performance of graphitic materials 
in any flow regime. 
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(3) Predictions with respect to integrated mass transfer and to 
shape change for basic geometric configurations, both in the stationary 
and in the spinning modes relative to the flow direction (flight path) 
can be obtained within the limits normally required for engineering 
accuracy. The multidimensional computer program (and with special 
care, the one dimensional monotonically shrinking sphere analysis) 
appear to be well suited for predicting the flight and mass transfer 
characteristics of re-entering objects using the theoretical-empirical 
expressions derived. However, specific flow environment and ex- 
perimental heat transfer distribution data must be factored into any 
computer analysis in order to adequately predict ground facility test 
performance. 
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APPENDIX A 

LOW DENSITY EFFECTS ON SPECIFICATION OF CHEMICAL KINETICS 
If one considers the oxidation reaction as first order one writes: 


and: 


m = k 
ox w 


( p ° 2 ) 


a) 


eff. 


If one considers eq. ( 1 ) in the limit of high Reynolds number: 
= .21 0 V 2 

05 00 

eff. 

2y 


( p ° 2 ) 


( 2 ) 


7 T C 


where y is a reaction probability. (3) 


However, if one attempts the use of eq. (1) in free molecular 
flow some modification is necessary in the definition of p^ , in this 


limit it is necessary to write for ^p^ ^ 


eff. 


K) 


eff. 


.21 P V 

<30 00 



- 

C 


If 

w 


2 

V 



00 


(4) 


and in the near free molecular regime: 



r 



II 

r 

O 

a 

P °2 


1 + . 5 Re 

00 

eff. 

eff. J 

T7* » 1 



F.M. 


(5) 


Using the results of Levinsky Ref, 19 for p in the low Reynolds number 

°2 

regime and interpolating between these results and eqs. (4) and (5) one 

obtains the results shown in Fig. 23 . 

The result for /p^ \ shown in Fig. 23 , can then be utilized 
' 2 'eff. 

in eq. ( 1 ) to obtain the mass transfer for reaction controlled oxidation. 
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APPENDIX B 

DERIVATION OF COMBINED 
THEORETICAL-EMPIRICAL EQUATIONS 


In order to graphically represent the combined theoretical-empirical 
mass loss performance of the materials studied in this program, the 
following expressions were employed: 


Let m^ 




-E/RT 

e 


represent rate controlled (1) 
mass loss 


Let m = C yj P /R 
D e rS 


represent diffusion con- (2) 
trolled mass loss 


Let m 



represent transition (3) 
from rate controlled 
to diffusion controlled 
mass loss 


Substitution of (1) and (2) into (3) yields, after some algebraic re -arrange - 



(4) 


from which the curves on Figures 32, 33, 52 and 53 were calculated. 

Further, if one lets 

f = fractional ratio of Stanton No. 

f = fractional ratio of effective pressure ratio 
2 

replacing m^ by f^ m^ and m^ by f^m^ in e< l» (3) one may derive 


f 2 k o 

r 2 E/RT A /f.k\“? 


L 2 E/KT + /t 2 K \ - R j 

1/2 

i V f i c / J 
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For the normalized mass loss rates a similar series of substitutions 


gave 


m 




k e 
o 


-E/RT 


.*» v-#/ • J 

for curves expressed in terms of the theory (Figures 35 and 54). 

The theory curves modified to correspond to data points acquired 
in their actual environments were determined from 

- E/RT 


2 -2 E/RT 

e 


i 


1/2 


(6) 





f k e 


J- >/ C Z “o\ 

2 -2 E/RT 1 

e 

1/2 

>B (f, C ) 




(7) 


This arrangement was used to relate the ATJ data more clearly with 
respect to the theory (Figure 34). 

In order to correlate data points with the theory curves (as repre- 
sented by eq'ns. (4) and (6), a slightly different approach was used in 
which it was necessary to determine a weighting factor for each test 
point. To establish this factor a tabulation of m was made for conditions 
under which the ratio of assumed particular values. For each 

material, the ratio 


m R 

rn 


k yR 

o v 


-E/RT 


B 


was evaluated and the appropriate range of temperatures was determined 
for each size specimen. The extent to which each mass rate (m^,m^) 
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contributed to the combined mass rate (m) for the corresponding condition 
was evaluated as the weighting factor for the datum point. See Appendix C. 


For a given set of nominal operating conditions, the factors f ^ and 
f^ were dependent upon the shock Reynolds number and hence were 
treated as constants to which the weighting factors were applied. 

In the preceding equations: 



ATJ Graphite 

Material "X" 

2 1/2 

k (lb. /sec. -ft -atm ) 

o 

4. 42 x 10 5 

2.09 x 10 3 

E (cal/mole) 

3 

44. 0 x 1 0 

3 

39. 1 x 1 0 

r nk , ,3/2 1/2 

C (lb. /sec. -ft atm ) 

-3 

6. 2 x 10 

8.5 x 10" 3 
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APPENDIX C 

DISCUSSION AND PRESENTATION OF SAMPLE CALCULATIONS 
REGARDING DATA REDUCTION 


To illustrate the process of reducing mass rate data in the experimental 
program, the results of two test runs are presented and carried through 
each of the calculations used. Where necessary, details of the specific 
procedures are expanded. 

Consider the results of runs 521 and 560 on 1/16” diameter ATJ 
graphite mass rate specimens (Table 15): 


Run 

No. 

Time 

sec. 

Length Chg. 
inches 

Mass Loss 
lb/ft. 

Mass Ra^e 
lb/ sec. ft. 

Ave« Surf. 
Temp. 
°K 

560 

250 

. 0750 

.6751 

— 

— 

521 

60 

.0166 

. 1494 

— 

— 


190 

.0584 

.5257 

2. 77 X 10" 3 

1515 


Difference in time and mass loss between the two runs at the same 

test conditions result in one mass rate data point at an average surface 

temperature (TJ of 1515°K. 

_ Am pA l 108 x, 0584 _ . 5257 

m A t 12 At 12 X A t At 

3 2 

m 1 = 2. 77 x 10” lb/sec. ft. (data pt. Figure 32, sym. # ) 
Other data, useful later in the computations are: 
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Reynolds number per foot = 520 ft, 

Reynolds number, Re § = = 2.7 

11 -3 

Effective radius, R^ = 2.46 x — x — = 6,41 x 10 ft. 

-3 

Model stagnation pressure, = 7. 70 x 10 atmos. 

The individual data points (m 1 ) were located on an Arrhenius plot 
(m* vs 1/T, Figure 32). Collectively, they were used to establish a 
reaction rate slope from which the specific reactivity (k^) an d activation 
energy (E) of the ATJ graphite were evaluated. The values obtained were 
expressed in a reaction rate equation, which when coupled with the theo- 
retical expression for diffusion controlled activity by the transition equation, 
provided a single equation (Appendix B) from which "theoretical" performance 
curves were evaluated. These were overplotted on the Arrhenius graph 
(Figure 32). 

The computed mass rates are for data acquired in regimes other 
than those of high Reynolds number, for whichthe original theory was 
derived. Hence experimental data have to be modified to correspond to 
the theory. The value of the modification factor (F^) depends on the 
degree to which oxidation is rate controlled or diffusion controlled* The 
plotted theory curves depict graphically a qualitative sense of the contri- 
bution of each process and indicate that the modifying factor required by 
either process must be weighted against that of the other process depending 
on their relative predominance for a given test point. The following dis- 
cussion constitutes the basis for obtaining a '’weight" applicable to a 
given data point. Every point was individually weighted. 
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1 


. 2 

m D 


1 


. 2 

m R 


1/2 (Transition equation, Sec. IIC, 
eq’n 30) 


Let N = m_ / m_ (o ^ N ^ c0 ). Then 
R D 


m = 

r l 

1 1 - 1 

'1/2 = 

1 

4 . 

1 


. 2 
[ m D 

M 2 * 2 

N m D j 



N 2 


1/2 


Divide both sides by m 


D 


m 


N 


II 

r , ii 

1/2 

,i 


1 + 2 


N +1 


N 2 




D 


The ratio gives o ^ 


1/2 


^ 1, which has the virtue of excluding 


explicit use of mass rates in the right hand member,, For any given 

value of N, the contributing mass rates are m and m = Nm . 

D R D 

The relative contribution of either process can be seen in the 
expression: 

1 


m 

~© 

rn 


D 


1 


_L ' 

N 2 


1/2 


If m is N times m (as postulated), it contributes — - toward the 
R D 

value of m when contributes unity, relative to it. Their total contri- 


bution is 1 + 


and their fractional contributions are 


hr 


for m R ; 


for m D ; 


1/N 


1 + 1/N 


1 + 1/N 


1 + N 
2 


N 


1 + N 
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Note that for N -» 00 , ratio m/m^ 
N -* 0, ratio m/m^ 
N = 1, ratio m/ 


1. (diff. control) 

N (rate control) 

\ yfl (equal control) 

u 


Numerical solutions follow: 



N = 


V -E/RT 

k V P e 




where C is the diffusion control level in lb/ft. 


3/2 1/2 

-sec -atm 


Then N 


N 


_ , -22150/1515) 

4. 42 x 10 (e ' 

-3 1/2 

6.2 x 10 (384) 

7.15 x 10 7 

19.6 (2. 705 x 10 6 ) 


N = 1.35 

Thus, the contributed fractions are 


for m R ; 


1 + N 


= .335 


f ° r ™d ; 


N 


1 + N 


= .645; where. 645 + . 355 = 1.00 


For the l/l6 u diameter mass rate models, Re g = 2.7; using this 
value one obtains a Stanton Number ratio (f^) of 0. 60 (Figure 22) for the 
diffusion controlled oxidation and an effective pressure ratio (f^) of 0.44 
(Figure 23) for rate controlled oxidation. 
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f i N 


- = F m ; Weighted Modifying Factor 


1 + N 1 + N 
0.60 (0. 645) + 0.44(0. 355) = 0. 543 
The originally calculated data point can now be modified, mf/F^ = m. 


thus, m = 


2. 77 x 10 
0. 543 


-3 


-3 2 

5.22 x 10 lb/ sec. ft. (data pt. , Figure 33) 


giving the value of the mass rate corresponding to that predicted by the 

original theory. The value of m 1 or of m may be normalized with respect 

to body radius and pressure and plotted against surface temperature (for 

this pt. , T = 1515°K). 
w 


Therefore, 


wS 


RB = Z.TTxlO - 3 / 6 ' 41 * 10 ' 3 ' 1/2 


7. 70 x 10 


-3 


-3 3/2 1/2 

= 2.54x10 lb/ft -sec-atm (data pt. Figure 34) 


Similarly, 




m 



B 


= 5.22 x 10 


-3 / 6.41 x 10 


-3 


1/2 


7. 70 x 10 


-3 


-3 3/2 1/2 

= 4.78x10 lb/ft -sec-atm (data pt. Figure 35) 
Thus, each datum point for comparison with the theory is established, 
calculated in the form pertinent to the type presentation for which it is to be 
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APPENDIX D 


DIFFUSION CONTROL LEVEL FOR MATERIAL "X" 

In the absence of available information regarding the chemical 
reactions likely to exist on the ablating surface of Material "X", a 
sound theoretical evaluation of the mass rate representative of 
diffusion controlled reactions cannot be obtained. The problem is 
probably further complicated by the presence of at least two reacting 
substances which may be combined with the available air to form 
various species. In view of the difficulties of attacking the problem 
in a strictly analytical manner, a method is presented in order to 
establish some criterion with which to compare the experimental data 
and to provide a conservative value for application to the systems 
analyses . 

Assuming the constituents to be uniformly and homogeneously 
distributed within the specimens, and postulating no preferential 
reactivity of oxygen with the constituents of Material "X", the mass 
fraction and number density relationships predict a normalized 
diffusion controlled mass rate of rn\J R^/P e = C-^=7.67 x 10 
Given n ^=n^ + n^ =1.0; where n^ = . 9875 and n^ = .0125 (Ref. 22) , 
and = n Q =1.0, one can weigh the contributions of each con- 


stituent to the normalized mass rate. 


n c m c + n u m u 


n c m c 


( catj ) 


UNCLASSIFIED 



then C = .9875 (12) + . 0125 (238) (6.2xl0 -3 ) 

-X. * 

1. 00 ( 12 ) 

= 7. 67 x 10' 3 

However, there is considerable evidence available from the 
experimental results which suggest that the uranium did indeed ex- 
hibit a greater affinity for the available oxygen. It was observed that 
when initially exposed to the hot air flow, there was a sizeable 
shower of brillant particles shed by the specimens, which gradually 
diminished with time. In addition, the motion pictures show local 
minute areas on the surface which become much more intensely 
bright than the overall surface, last for a short time and then dis- 
appear. These phenomena recall the well-known behavior of uranium 
and uranium carbides having small specific areas to ignite spontan- 
eously in air at very moderate temperatures with a very significant 
generation of heat (Ref. 23). Close examination of the specimen 
surfaces of Material^X" tends to confirm that such a process does 
occur. The fissured and porous appearance of the specimen surfaces 
after exposure suggest that indeed, the oxygen combines preferentially 
with the uranium. 

Without well defined data to establish the degree to which a select- 
ivity phenomenon might operate, it was postulated from an analysis 
of the data in Table 4 that a conservative value might be a rate 1-1/2 
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times the parity loss value of uranium (see Section III E # 5). Cal- 
culated on this basis, 


'X 


m r 1 1.5 n TT m TT / C ^ 

c C U u ( ATJ ) 


n c m c 


one finds the diffusion controlled mass rate parameter C ~ 8. 5x 10 . 

X 

This value was then combined with the empirically evaluated reaction 

rate expression using the transition equation (Section II C, equation 30), 

and empirical-theory curves were calculated and plotted with the 

mass rate experimental data (Figure 54). Examination of the graph 

clearly indicates that the experimental data define no trend toward a 

leveling off of mass loss over the range explored. Rather, the data 

points suggest a continuance of the exponential rise of mass rate 

-3 

parameter up to 10.5 x 10 

The results of the experiment indicate that selective reactions do 
occur and, indeed that oxygen has, perhaps, an even greater affinity 
for the uranium constituent of Material "X" than is postulated. 
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APPENDIX E 

COMPARISON OF EXACT AND NUMERICAL 
HEAT CONDUCTION SOLUTIONS FOR SPHERES 


The transient heat conduction equation for a sphere heated by a 
constant rate is 


„ ST _ 1 S / 2 ST \ 

0C p St 2 S R ( KR BR/ 

R 


(i) 


and the following boundary and initial conditions applied: 

T(R,o) = T e 

ST °w 

TR (a>t) = IT 


( 2 ) 


ST 

SR 


(o.t) = 0 


The exact solution to this problem is given in Reference 25 as: 


T - T 


3Q t Q (5R 2 - 3a 2 ) 
w + w 


)0C a 
P 


lOKa 


a i t 


a ' l 


R 


2Q w a 2 ^ sin (-^) 

r 2 • >• 

n=l C sin £ 


KR 


(3) 


where is the uniform heating rate 

a is the radius of the sphere 

£ (n = 1, 2, 3...) are the positive real roots of the 
transcendental equation, tan £ = £ 

is the thermal diffusitivity = ( ) 

DC p 
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For R = a, equation (3) reduces to 


2 

C at 

n 


T - T 


3Q t 
w 

DCp a 


Q a 

W 

5K 


2Q a 
w 

IT“ 


00 . yj „ a C 

sin C n £ a 

2 , ~2 . 
n=l C n smC n 


( 4 ) 


Figure 103 illustrated the comparison of the surface temperature as a 
function of time obtained from the IBM 7094 computer solution and the 
exact solution of equation (4). At 3000 R, the percent error of the 
numerical technique is less than 2% for the case considered. Further 
improvement in accuracy is possible for this comparison by decreasing 
the time interval between computations. 

The transient heat conduction equation for a sphere analyzed muli- 
dimensionally is 


pc 


5T 1 a / 2 9T\ 

at 2 dR V SR/ 
R 


R sin B 


b 

be 




+ 


1 

2 . 2 . 

R sm 0 



(5) 


In Reference 26, a three dimensional solution was obtained for the case 

o 

of a sphere rotating about an axis 90 to the flow direction. If the spin 
rate is set to zero, the solution may be applied to a sphere whose axis 
of rotation is aligned to the flow direction and the following boundary 
applied 

K (a, 6, (p) = (1 + sin 0cos <p) 

oR ^ 
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where the product, sin 0 cos <p = constant, locates a plane perpendicular 

to the flow direction (x - direction) as illustrated in the following sketch. 

z 



Since the flow is axisymmetr ic, the solution is obtained for a 
specific heat flux distribution in the x-y plane, Figure 104, where 

Q 


K I< a ’ !• * 


w o 


( 1 + cos (p) 


(7) 


The surface temperature T is given as the sum of four parts: 

T = T_ + T. + T + T where 
E 1 l , 5 

T^, is the initial temperature, is the solution obtained for a constant 
heat flux, — as listed in Equation (3) and 


Q R 

T 2 + T 3 sin9cos ‘p 


R 2a 

a + _2 — 

R n=l 


2 


c at 

n 

R 2 




/ . «, R *> R R \ 

( sinf — - C — cos C — ) 
\ n a n a n a / 


C cosC 

n n 
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Figure 104 illustrates the specific surface heat flux distribution for which an 
exact analytic solution exists. Figure 105 compares this exact solution with the 
numerical solution obtained from the IBM 7094 computer program at three 
surface locations (R = a) in order to check the accuracy of the numerical 
solution. The numerical solution deviates by less than 3/4 of 1% at 4000°R. 
Note that this chosen heat flux distribution does not pertain to the systems 
analyses per se. 
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APPENDIX F 

DETERMINATION OF AVERAGE HEAT AND MASS 
TRANSFER FOR SPHERES AND CYLINDERS 


In order to analyze the sphere one -dimensionally, a suitably 
averaged heat and mass transfer rate must be determined. Assuming 
a Newtonian pressure distribution on a hemisphere, the heat and mass 
transfer at a given body angle, becomes 


Q 


= Q 


local 


w cos 0 
o 


( 1 ) 


i = m n 

w „ w cos 0 

local o 


As a first approximation, one may base the average heat and mass 
transfer on a surface -area basis: 


Q 


= /^local dA 


hemisphere 


or in spherical coordinates 


w 


hemisphere 


/dA 

A 


(2) 



( 3 ) 


hemisphere = . 5Q 


( 4 ) 
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If one assumes the aft heating (and mass transfer) to be zero, then for 
the entire sphere, the average heat transfer would be 

“Q = .25 Q (5) 

w , w 

sphere o 

For the systems study an average value between . 5 and .25 of 

the stagnation heating rate (and mass transfer), viz. , . 375 was assumed. 

It should be noted that a more exact calculation for the average heat 

transfer to a non-ablating hemisphere was found in Reference 2 7 to be 

0. 436 Q . The effects of varying the averaging factor between . 5 and 
o 

.2 were investigated in the systems study (Figure 82). For the cylindrical 

rod an average value equal to the stagnation value was imposed over the 
£/hx. 

entire front iaot < i 
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TABLE 1 

SUMMARY OF PRESSURE PROBE MEASUREMENTS 


Position * 

3 

Stagnation Pressure, Atmospheres x 10 

High En 

thalov 

Low Enthalpy 

5. 0" Station 

1.2" Station 

5. 0" Station 

0.563 

8. 25 

56. 0 

7.00 

1.000 

7.98 



1.500 

7.70 




* Distance from Nozzle Exit 

Note: Each tabulated value represents the average of the several test 
exposures made at that position. Maximum variation in results 
was + 3. 5% 
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MASS LOSS DATA 
1 /4" Dia. Cylinders 
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(1) See Fig. 12, Item 17 

(2) F osition: 0. 5625" from nozzle exit 

(3) See Fig.l 3, (c) 

(4) For two color surface temperature data, see Table 11 
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UNCLASSIFIED 


(1) See Fig.12, Item 21 

(2) Fosition: 0. 5625 n from nozzle exit 

(3) See Fig. 13, (b) 

(4) See Fig. 13, (a) 

(5) For two color surface temperature data, see Table 12 
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(1) See Fig, 12, Item 21 

(2) Fosition: 0, 5625" from nozzle exit 

(3) Fosition: 1,50" from nozzle exit 

(4) See Fig, 1 3 , (d) 

( 5 ) For two color surface temperature data, see Table 15 
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(1) See Fig. 12, Items 13, 19 

(2) Position: 0, 800 M from nozzle exit 

(3) For two color surface temperature data, see Table 16 



TABLE 3 
MASS LOSS DATA 
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(1) See Fig. 12, Items 18, 19 

(2) Position: 0.800" from nozzle exit 
(#) Optical pyrometer 

(3) For two color surface temperature data, see Table 21 
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TABLE 4 

COMPARISON: MASS LOSS DATA 
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ONE INC H CYLINDERS WITH 1/4" INSERTS 


Length Los 

s Basis 

Insert Mass Loss Basis 

Ratio * 

length loss 

mass loss 

mass loss 

mass loss 
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lb/ ft. 2 

grams 

. lb/ ft. 2 
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. 1902 

1. 230 
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UNCLASSIFIED 


(1) See Fig. 12 Item 13 (4) Position:!. 85(7* from nozzle exit (7) Air Mass Flow .0014 lb/sec 

(2) See Fig. 12 Item 14 (5) Position: 2. 225" from nozzle exit (*) Optical pyrometer 

(3) Position: L. 500" from nozzle exit (6) Position: 0. 800^' from nozzle exit 








INTEGRATED MASS LOSS DATA 



m in in mm 
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oomooomooc~ooo^HC-oo 

men— ^MnrOHNHCOTfNro^o^ 


POOOin^hHOrnGO^NrOO(>NsO 
COCOCO— (COCOCO— I— trncOCOCOOvOvD 
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(1 ) See Fig.lZjtem 13 (4) Position: 1 . 850" from nozzle exit (7) Air Mass Flow .0014 lb/ 

(2) See Fig.l2,Item 14 (5) Position: 2. 225" from nozzle exit (*) Optical pyrometer 

(3) Position: 1.500" from nozzle exit (6) Position: 0. 800" from nozzle exit 
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TABLE 8 

SUMMARY OF TOTAL CALORIMETER 
ENTHALPY MEASUREMENTS 


E k ,KW 

iti'lb/sec. 

P s,psia 

fcs/ RT o 

h s> Btu /lb. 

238 

.0 

)15 

22.2 

351 

11880 

238 



20.8 

327 

11070 

230 



20.6 

316 

10700 

283 


f 

19.5 

332 

11230 

76.7 

.0( 

114 

16.0 

157 

5310 

76.7 



15.9 

166 

5610 

76.7 



16.0 

162 

5480 

87.5 

> 

f 

16. 3 

162 

5480 
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UNCLASSIFIED 









SUMMARY OF COLD WALL HEAT TRANSFER MEASUREMENTS 



PC-LASSTO 





UNCLASSIFIED 


(1) Configuration: Gives model dia. followed by copper insert dia< 

(2) Station: Distance to stagnation point, from nozzle exit 

(3) Copper cylinder with .003" RTV coating 








SUMMARY OF COLD WALL HEAT TRANSFER MEASUREMENTS 
SPHERICAL SPECIMENS 





UNCLASSIFIED 


(1 ) Configuration: Gives model dia. followed by copper insert dia. (2) Station: Distance to stag, point, from nozzle exit 

(3) Decayed model position increased by amount ablated @ stag, point 

(4) These values corrected for relatively oversize insert-see text 
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TABLE 11 

SUMMARY OF MASS RATE DATA 


Specimen Dia. 

1/4" 

Notes - 

Specimen Mat’l 

ATJ 

(1) 

Nozzle Exit Dia. 

5" 

(2) 

Stag’n Press. , atm 

8. 25 x 10~ 3 


Enthalpy Ratio 

350 

0 

Reynolds No. , Re g 

11 



Unmodified data 
Data modified for 
Reynolds No. Effects 
Symbol on graphs 
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TABLE 1 1 

SUMMARY OF MASS RATE DATA 


Specimen Dia. 

1/4" 

Notes - 


Specimen Mat'l 

ATJ 

(1) 

Unmodified data 

Nozzle Exit Dia. 

5" 

(2) 

Data modified for 

Stag'n Press. , atm 

8. 25 x 10- 3 


Reynolds No. Effects 

Enthalpy Ratio 

350 

0 

Symbol on graphs 

Reynolds No., Re 

11 
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TABLE 12 


SUMMARY OF MASS RATE DATA 


Specimen Dia. 

1/4” 

Note* - 


Specimen Mat'l 

ATJ 

(1) 

Unmodified data 

Nozzle Exit Dia. 

1.2" 

(2) 

Data modified for 

Stag'n Press. , atm 

5. 60 x 10~ 2 

© 

Reynolds No. Effects 

Enthalpy Ratio 

350 

Symbol on graphs 

Reynolds No. , Re 

85 
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TABLE 12 

SUMMARY OF MASS RATE DATA 


Specimen Dia. 

1/4” 

Notes - 


Specimen Mat 1 ! 

ATJ 

(1) 

Unmodified data 

Nozzle Exit Dia, 

1.2” 

(2) 

Data modified for 

Stag'n Press, , atm 

5. 60 x 10~ 2 


Reynolds No. Effects 

Enthalpy Ratio 

350 

© 

Symbol on graphs 

Reynold^ No. , Re 

85 
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TABLE 12 


SUMMARY- 0F MASS RATE DATA 


Specimen Dia. 

1/4 M 

Notes - 


Specimen Mat'l 

ATJ 

(1) 

Unmodified data 

Nozzle Exit Dia. 

1. 2" 

(2) 

Data modified for 

Stag'n Press. , atm 

5. 60 x 10”^ 


Reynolds No. Effects 

Enthalpy Ratio 

350 

© 

Symbol on graphs 

Reynolds No. , Re 

s 

85 
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TABLE 13 


SUMMARY OF MASS RATE DATA 


Specimen Dia, 

1" 

Notes - 


Specimen Mat'l 

ATJ 

(1) 

Unmodified data 

Nozzle Exit Dia. 

5" 

(2) 

Data modified for 

Stag'n Press. , atm 

9.7 x ICT 3 


Reynolds No. Effects 

Enthalpy Ratio 

400 

A 

Symbol on graphs 

Reynolds No., Re g 

44 
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TABLE 1 3 

SUMMARY OF MASS RATE DATA 


Specimen Dia. 

1" 

Note* - 


Specimen M&t’l 

ATJ 

(1) 

Unmodified data 

Nozzle Exit Dia. 

5" 

(2) 

Data modified for 

Stag'n Press. , atm 

9. 7 x 10’ 3 

A 

Reynolds No. Effects 

Enthalpy Ratio 

400 

Symbol on graphs 

Reynolds No. , Re g 

44 
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TABLE 13 

SUMMARY OF MASS RATE DATA 


Specimen Dia. 


Notea - 


Specimen Mat'l 

ATJ 

(1) 

Unmodified data 

Nozzle Exit Dia. 

5" 

(2) 

Data modified for 

Stag'n Press. , atm 

9.7 x 10’ 3 


Reynolds No. Effects 

Enthalpy Ratio 

400 

A 

Symbol on graphs 

Reynolds No. , Re g 

44 
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TABLE 1 3 

SUMMARY OF MASS RATE DATA 


Specimen Dia, 

1" 

Notes - 


Specimen Mat'l 

ATJ 

(1) 

Unmodified data 

Nozzle Exit Dia. 

5 M 

(2) 

Data modified for 

Stag'n Press. , atm 

9.7 X 1CT 3 


Reynolds No. Effects 

Enthalpy Ratio 

400 

A 

Symbol on graphs 

Reynolds No., Re g 

44 
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TABLE 14 

SUMMARY OF MASS RATE DATA 


Specimen Dia, 

1" 

Notes - 

Specimen Mat'l 

ATJ 

(1) Unmodified data 

Nozzle Exit Dia. 

1.2" 

(2) Data modified for 

Stag'n Press. , atm 

4.43 x 1CT 2 

Reynolds No. Effects 

Enthalpy Ratio 

400 

^ Symbol on graphs 

Reynolds No. , Re g 

342 
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TABLE 15 


SUMMARY OF MASS RATE DATA 


Specimen Dia. 

1/16" 

Notes - 


Specimen Mat ( l 

ATJ 

(1) 

Unmodified data 

Nozzle Exit Dia. 

5" 

(2) 

Data modified for 

Stag'n Press* , atm 

7.70 x 10~ 3 


Reynolds No. Effects 

Enthalpy Ratio 

Reynolds No. , Re^ 

350 

2.7 

• 

Symbol on graphs 




Mass Loss 
Lb/ ft. 2 

Avg. 

Mass Rate it 10 3 

• XE 

m \l p« x 

10 3 

Run 

No*. 

Time 

Sec. 

°K 

Temp. 

lb/ft 2 - 

sec 

lb/ft y2 -*ec-Atm 1/2 

(1) 

(2) 

(1) 

(2) 

560 

n 

. 6751 






521 

KSp 

. 1494 








. 5257 

1515 

2. 77 

5. 21 

2. 54 

4. 78 

560 

1 

. 6751 






520 

HI: 

.3168 







130 

. 3583 

1 

1520 

2. 75 

5. 15 

2. 52 

4. 73 

560 

Em 

. 6751 






519 

BM 

. 2988 







130 

. 3763 

1540 

2.90 

5. 38 

2. 66 

4. 94 

560 

250 

. 6751 






518 

100 

. 2556 







150 

.4195 

1545 

2. 80 

5. 17 

2. 56 

4. 73 

560 

250 

. 6751 






517 

100 

. 2520 







150 

.4231 

1560 

2. 82 

5. 15 

2. 58 

4. 72 

560 

250 

. 6751 






516 

80 

.1863 







170 

.4888 

1545 

2.87 

5. 30 

2. 62 

4. 84 

560 

250 

. 6751 






515 

80 

. 2007 







170 

.4744 

1520 

2. 79 

5. 23 

2. 55 

4. 79 

560 


. 6751 






514 

60 

. 1620 







190 

. 5131 

: 1500 

2. 70 

5. 13 

2.47 

4. 70 

560 


. 6751 






513 


.1107 







210 

. 5644 

1485 

2. 69 

5. 17 

2.46 

4. 72 
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TABLE 15 

SUMMARY OF MASS RATE DATA 


Spe cimen Dia. 1 / 16 " 

Specimen Mat*l ATJ 

Nozzle Exit Dia. 5 " 

Stag'n Press., atm 7 # 70 x 10“ 3 
Enthalpy Ratio 350 

Reynolds No. , Re^ 2.7 


Notes - 

( 1 ) Unmodified data 

( 2 ) Data modified for 
Reynolds No. Effects 

• Symbol on graphs 
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Specimen Dia. 
Specimen Mat'l 
Nozzle Exit Dia. 
Stag'n Press. f atm 
Enthalpy Ratio 
Reynolds No. , Re 

s 


TABLE 15 

SUMMARY OF MASS RATE DATA 
1/16" Notes - 

ATJ (1) 

5” „ (2) 


7.70 x 

350 

2.7 


10 


-3 


Unmodified data 
Data modified for 
Reynolds No. Effects 
Symbollon graphs 
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TABLE 15 

SUMMARY OF MASS RATE DATA 


Spe cimen Dia. 1 / 1 6 n 

Specimen Mat*l ATJ 

Nozzle Exit Dia* 5" 

Stag'n Press. , atm 7 # 70 x 10~ 3 

Enthalpy Ratio 35 0 

Reynolds No., Re 2.7 

s 


Notes * 

(1) Unmodified data 

(2) Data modified for 
Reynolds No. Effects 

% Symbol on graphs 
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TABLE 15 

SUMMARY OF MASS RATE DATA 

Specimen Dia. 1/16" Notes - 

Specimen Mat' 1 ATJ ( 1 ) Unmodified data 

Nozzle Exit Dia. 5" (2) Data modified for 

Stag'n Press., atm 7. 70 x 10' J Reynolds No. Effects 

Enthalpy Ratio 350 • Symbol on graphs 

Reynolds No. , Re 2. 7 

s 


Run 

Nos. 

Time 

Sec. 

Mass Loss 
Lb/ ft. 2 

Avg. 

°K 

Temp. 

520 

120 

. 3168 


513 

40 

.1107 



80 

. 2061 

1435 

519 

120 

. 2988 


521 

60 

. 1494 



60 

. 1494 

1480 

519 

120 

. 2988 


518 

100 

. 2556 



20 

. 0432 

1510 

519 

120 

. 2988 


517 

100 

. 2520 



20 

. 0468 

1475 

519 

120 

. 2988 


516 

80 

. 1863 



40 

.1125 

1510 

519 

120 

. 2988 


515 

80 

. 2007 



40 

. 0979 

1485 

519 

120 

.2988 


514 

60 

. 1620 



60 

. 1368 

1465 

519 

120 

. 2988 


513 

40 

.1107 



80 

. 1881 

1450 

518 

100 

. 2556 


521 

60 

. 1494 



40 

. 1062 

1475 
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TABLE 15 

SUMMARY OF MASS RATE DATA 


Specimen Dia, 

1/16’' 

Notes - 


Specimen Mat'l 

ATJ 

(1) 

Unmodified data 

Nozzle Exit Dia. 

5" ^ 

(2) 

Data modified for 

Stag'n Press. , atm 

7. 70 x 10 


Reynolds No. Effects 

Enthalpy Ratio 

350 

• 

Symbol on graphs 

Reynolds No. , Re 

s 

2.7 



Run 

Nos. 

Time 

Sec. 

Mass Loss 
Lb /ft. 2 

Avg. 

°K 

Mass Rate k 10^ 
lb/ft 2 - sec 

* TE 3 

m \l Pe x 10 

lb/ft y2 -sec-Atm 1/2 

Temp. 

(1) 

vi 

(1) 

(2) 

518 

100 

. 2556 


| 




516 

80 

. 1863 







20 

. 0693 

1505 

3.46 

6. 55 

3. 16 

5. 99 

518 

100 

. 2556 






515 

80 

. 2007 







20 

. 0549 

1480 

2. 74 

5. 20 

2. 51 

4. 84 

518 

100 

. 2556 






514 

60 

. 1620 







40 

. 0936 

1460 

2. 34 

4. 58 

2. 14 

4.17 

518 

100 

. 2556 






513 

40 

.1107 







60 

.1449 

1445 

2.41 

4. 75 

2. 21 

4. 36 

517 

100 

. 2520 






521 

60 

. 1494 







40 

.1026 

1495 

2. 56 

4. 89 

2. 34 

4. 47 

517 

100 

. 2520 






516 

80 

. 1494 







20 

. 0657 

1525 

3. 29 

6. 15 

3. 01 

5. 63 

517 

100 

. 2520 






515 

80 

. 2007 







20 

.0513 

1500 

2. 56 

4. 85 

2. 34 

4.45 

517 

100 

.2520 






514 

60 

. 1620 







.40 

. 0900 

1480 

2. 24 

4. 26 

2. 05 

3.95 

517 

100 

. 2520 






513 

40 

. 1107 







60 

. 1413 

1465 

2. 35 

4. 59 

2. 15 

4. 19 


173 


UNCLASSIFIED 


i-uiwUssiliio::- 


• • • • • 

• • • • • • • 

• • • • • 

• •••••• 


TABLE 15 

SUMMARY OF MASS RATE DATA 


Specimen Dia. 

1/16" 

Notea - 


Specimen Mat'l 

ATJ 

(1) 

Unmodified data 

Nozzle Exit Dia. 

5" 

(2) 

Data modified for 

Stag'n Press. , atm 

7. 70 x 10' 3 


Reynolds No. Effects 

Enthalpy Ratio 

Reynolds No. , Re 

s 

250 
2: 7 

• 

Symbol on graphs 


Run 

Noa. 

Time 

Sec. 

Mass Loss 
Lb/ft. 2 

>o 

< 1 

Mass Rate k 10 3 
lb/ft 2 - aec 

* rs 3 

m \| Pe x i<r 

lb/ft^ 2 -sec-Atm 1/2 

Temp. 

(1) 

(2) 

(1) 

(2) 

516 

80 

. 1863 






521 

60 

. 1494 







20 

1 . 0369 

1470 

1. 85 

3. 60 

1.. 69 

3. 28 

516 

80 

. 1863 






513 

40 

.1107 







40 

. 0756 

1440 

1. 89 

3. 74 

1. 73 

3.43 

515 

80 

. 2007 






521 

60 

. 1494 







20 

. 0513 

1475 

2. 56 

5. 00 

2. 34 

4. 53 

515 

80 

. 2007 






514 

60 

. 1620 







20 

. 0387 

1460 

1. 93 

3. 77 

1. 77 

3.46 

515 

80 

. 2007 






513 

40 

. 1107 







40 

. 0900 

1445 

2. 24 

4.42 

2. 05 

4. 04 

521 

60 

. 1494 






513 

40 

. 1107 







20 

.0387 

1445 

1. 94 

3. 82 

1. 78 

3. 51 

514 

60 

.1620 






513 

40 

.1107 







20 

.0513 

1445 

2. 56 

5. 05 

2. 34 

4. 62 

521 

60 

; 1494 

1400 

2. 49 

5. 09 

2. 28 

4. 65 

516 

80 

.1863 

1460 

2. 33 

4. 55 

2. 13 

4. 17 

520 

120 

.3168 

1410 

2. 64 

5. 35 

2.42 

4. 90 
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Specimen Dia. 
Specimen Mat'l 
Nozzle Exit Dia. 
Stag’n Press. , atm 
Enthalpy Ratio 
Reynolds No. , Re 

s 


TABLE 16 

SUMMARY OF MASS RATE DATA 


1 " 

ATJ 

5" 

7. 00 x 10 
160 
51 


Notes - 

(1) Unmodified data 

(2) Data modified for 
Reynolds No. Effects 

^ Symbol pn graphs 


Run 

No«. 

Time 

Sec. 

Mass Loss 
Lb/ ft. 2 

Avg. 

°K 

Mass Rate k 10^ 
lb/ft ^ - sec 

• .in' 3 

m M Pe x 10 

lb/ft^-sec-Atm^ 2 

Temp. 

cn 


(1) 

(2) 

631 

1200 

1.2195 

1280 

1. 02 

1. 042 

3. 90 

3. 98 

631 

1200 

1. 2195 






634 

240 

. 2232 







960 

.9963 

1 

1270 

1. 04 

1. 09 

3. 975 

4. 165 

631 

1200 

1. 2195 






642 

310 

. 2305 







890 

. 9890 

1280 

1.11 

1. 134 

4. 24 

4. 34 

631 

1200 

1. 2195 






635 

360 

. 2979 







840 

. 9216 

1285 

1. 10 

1. 123 

4. 21 

4. 29 

631 

1200 

1. 2195 






633 

463 

. 4626 







737 

. 7569 

1260 

1. 03 

1. 090 

3. 94 

4. 165 

631 

1200 

1. 2195 






630 

600 

.4707 







600 

. 7488 

1285 

1. 25 

1. 278 

4. 77 

4. 89 

631 

1200 

1.2195 






638 

720 

.6138 







480 

. 6057 

1275 

1. 26 

1. 307 

4. 82 

5. 00 

631 

1200 

1. 2195 






643 

772 

. 6462 







428 

. . 5733 

1300 

1. 34 

1. 348 

5. 12 

5. 15 

631 

1200 

1/2195 






644 

801 

. 7137 







399 

. 5058 

1295 

1. 265 

1. 287 

4. 84 

4.92 

632 

960 

. 6930 

1260 

0. 722 

. 765 

2. 76 

2. 92 
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TABLE 16 

SUMMARY OF MASS RATE DATA 


Specimen Dia. 

1” 

Specimen Mat*l 

ATJ 

Nozzle Exit Dia. 

5" 

Stag'n Press. , atm 

7. 00 x 10 

Enthalpy Ratio 

160 

Reynolds No. , Re g 

51 


Notes - 

(1) Unmodified data 

(2) Data modified for 
Reynolds No. Effects 

^7 Symbol on graphs 
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TABLE 16 

SUMMARY OF MASS RATE DATA 


Specimen Dia. 

1" 

Notes - 


Specimen Mat'l 

ATJ 

(1) 

Unmodified data 

Nozzle Exit Dia. 

5M * 

(2) 

Data modified for 

Stag*n Press. f atm 

7. 00 x 10’° 


Reynolds No. Effects 

Enthalpy Ratio 

160 

V 

Symbol on graphs 

Reynolds No. , Re g 

51 




Run 

Noe. 

Time 

Sec. 

Mass Loss 
Lb/ ft. 2 

Avg. 

°K 

Temp. 

Mass Rate k 10 3 
lb/ft 2 - sec 

• T* 3 

m \| Pe x 10 J 

lb/ft' V2 -8ec-Atm 1/2 

a) 


(1) 

... 2) 

644 

801 

. 7137 






642 

310 

. 2305 







491 

.4832 

1255 

0. 985 

| 1.048 

3. 76 

4. 00 

644 

801 

: 7137 






635 

360 

. 2979 







441 

.4158 

1 340 

0. 944 

. 895 

3. 61 

3.42 

644 

801 

. 7137 






633 

463 

.4626 







338 

. 2511 

1 315 

0. 744 

. 733 

2. 85 

2. 91 

644 

801 

. 7137 






641 

480 

.4104 







321 

. 3033 

1270 

0. 942 

. 988 

3. 60 

3. 775 

644 

801 

. 7137 






639 

600 

.4707 







201 

. 24 30 

1340 

1. 210 

1. 147 

4. 62 

4. 38 

644 

801 

. 7137 






638 

720 

. 6138 







81 

.0999 

1330 

1. 233 

1. 193 

4. 71 

4. 56 

643 

772 

. 6462 

1 





634 

240 

. 2232 







532 

.4230 

1280 

. 795 

. 820 

3. 04 

3. 135 

643 

772 

.6462 






642 

310 

. 2305 







462 

.4157 

1285 

. 900 

. 919 

3.44 

3. 51 

643 

772 

. 6462 






635 

360 

.2979 







412 

. 3483 

1295 

. 845 

. 854 

3. 23 

3. 265 
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TABLE 16 

SUMMARY OF MASS RATE DATA 


Specimen Dia. 

1” 

Notes - 

Specimen Mat’l 

ATJ 

(1) 

Nozzle Exit Dia. 

5" 

(2) 

Stag'n Press. # atm 

7. 00 x 10’ 3 


Enthalpy Ratio 

160 

V 

Reynolds No. , Re g 

51 



Unmodified data 
Data modified for 
Reynolds No. Effects 

Symbol on graphs 
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TABLE 1 6 

SUMMARY OF MASS RATE DATA 


Specimen Dia. 

1 ” 

Notea - 


Specimen Mat'l 

ATJ 

(1) 

Unmodified data 

Nozzle Exit Dia. 

5" 

(2) 

Data modified for 

Stag'n Press. , atm 

7. 00 x 10 “ 3 


Reynolds No. Effects 

Enthalpy Ratio 

160 


Symbol on graphs 

Reynolds No. , Re, g 

51 




Run 

Noa. 

Time 

Sec. 

Mass Loss 
Lb/ ft. 2 

Avg. 

°K 

Mass Rate k 10 ^ 
lb/ft 2 - aec 

a . 3 

m \l Pe * 10 ° 

lb/ft y2 -aec-Atm 1/2 

Temp, 

TT7 

m 

( 1 ) 

( 2 ) 

638 

720 

. 6138 






639 

600 

.4707 







120 

. 1431 

1300 

1 . 190 

1. 197 

4. 55 

4. 57 

642 

310 

. 2305 

j 

1250 

0 . 745 

. 798 

2. 85 

3. 05 

639 

600 

.4707 

1245 

0. 785 

. 849 

3. 00 

3. 24 

639 

600 

.4707 






634 

240 

. 2232 







360 

. 2475 

1295 

0 . 688 

. 695 

2. 63 

2 , 66 

639 

600 

. 4707 






642 

310 

. 2305 ! 







290 

. 2402 

1225 

0.829 

. 917 

3. 17 

3. 51 

639 

600 

.4707 






635 

360 

. 2979 







240 

. 1728 

1320 

0. 720 

. 703 

2. 75 

2 . 69 

639 

600 

. 4707 






641 

480 

.4104 







120 

. 0603 

1240 

0. 504 

. 549 

1. 927 

2 . 10 

641 

480 

.4104 

1255 

0. 855 

. 911 

3. 27 

3.48 

641 

480 

.4104 






634 

240 

. 2232 







240 

a 1872 

1270 

0. 780 

. 817 

2 . 98 

3. 125 

641 

480 

.4104 






642 

3*0 

. 2305 







170 

.1799 

1275 

1. 068 

1.108 

4. 08 

4. 23 

641 

480 

.4104 






635 

360 

. 2979 







120 

i 

.1125 

1285 

0. 938 

. 958 

3. 585 

3. 66 
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TABLE 16 

SUMMARY OF MASS RATE DATA 


Specimen Dia. 

1 " 

Note. - 

Specimen Mat 1 ! 

ATJ 

(1) 

Nozzle Exit Dia. 

5 ” * 

(2) 

Stag'n Press. , atm 

7 . 00 x 10"- 5 


Enthalpy Ratio 

160 

V 

Reynolds No. , Re^ 

51 



Unmodified data 
Data modified for 
Reynolds No. Effects 
Symbol on graphs 
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UNCLASSIFIED 



SUMMARY OF MASS RATE DATA 


Specimen Dia. 

1/4” 

Notes - 


Specimen Mat'l 

Mat’l. X 

(1) 

Unmodified data 

Nozzle Exit Dia. 

5" 

(2) 

Data modified for 

Stag'n Press. , atm 

8. 25 x 10‘ 3 


Reynolds No. Effects 

Enthalpy Ratio 

350 

B 

Symbol on graphs 

Reynolds No. , Re g 

11 




Run 

Nos. 

Time 

Sec. 

Mass Loss 
Lb/ ft. 2 

Avg. 

°K 

3 

Mass Rate k 10 
lb/ft 2 - sec 

. v flE 3 

m \| Pe x 10 

lb/ft^-sec-Atm 1 ^ 2 

Temp. 

rn 


(i) 

(2) 

690 

300 

. 3705 

1650 

1. 235 

1.995 

2. 19 

3. 54 

690 

300 

. 3705 






688 

240 

. 2500 







60 

. 1205 

1765 

2.008 

3; 24 

3. 56 

5. 74 

690 

300 

. 3705 






685 

180 

. 1853 







120 

. 1852 

1700 

1. 54 

2.485 

2. 73 

4. 40 

690 

I 

300 

. 3705 






687 

175 

. 1676 







125 

. 2029 

1705 

1. 63 

2. 63 

2. 88 

4. 65 

690 

300 

. 3705 






689 

120 

.1218 







180 

. 2487 

1775 

1. 38 

2. 23 

2.44 

3. 94 

690 

300 

. 3705 






686 

60 

. 0502 







240 

. 3203 

1725 

1. 34 

2. 16 

2. 37 

3. 82 

688 

240 

. 2500 

1710 

1. 04 

1.688 

1. 84 

2. 97 

688 

240 

. 2500 






685 

180 

. 1853 







60 

. 0647 

1690 

1. 08 

1. 74 

1.91 

3. 08 

688 

240 

. . 2500 






687 

175 

.1676 







65 

. 0824 

1695 

1. 29 

2. 08 

2. 28 

3. 68 

688 

240 

. 2500 






689 

120 

.1218 







120 

. 1282 

1765 

1. 07 

1. 726 

1. 89 

3. 05 
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Specimen Dia. 
Specimen Mat'l 
Nozzle Exit Dia. 
Stag'n Press. , atm 
Enthalpy Ratio 
Reynolds No. , Re 


TABLE 17 

SUMMARY OF MASS RATE DATA 
1/4" Notes - 

Mat’l. X (1) 

5” (2) 

8. 25 x 10~ 3 
350 


0 


11 


Unmodified data 
Data modified for 
Reynolds No. Effects 

Symbol on graphs 













Specimen Dia. 
Specimen Mat'l 
Nozzle Exit Dia. 
Stag’n Press. , atm 
Enthalpy Ratio 
Reynolds No. , Re 



TABLE 17 

SUMMARY OF MASS RATE DATA 

1 / 4" Notes - 


Mat'l. X 
5 " 

8. 25 x 10" 3 
350 

11 


(1) Unmodified data 

(2) Data modified for 
Reynolds No. Effects 

0 Symbol on graphs 


Run 

Nos. 

Time 

Sec. 

Mass Loss 
Lb/ ft. 2 

w> & 

> 0 

< 

Mass Rate k 10^ 
lb/ft 2 - sec 

. J Kb 3 

m\J p e x 10 J 

lb/ft^-sec-Atm^ 

Temp. 

(I) 

(2) 

(i) 

(2) 

695 

360 

. 3970 






688 

240 

. 2500 







120 

. 1470 

1790 

1. 23 

1.985 

2. 18 

3. 52 

695 

360 

. 3970 






685 

180 

. 1853 







180 

.2117 

1725 

1.175 

1. 90 

2. 08 

3. 36 

695 

360 

. 3970 






687 

175 

. 1676 







185 

. 2294 

1730 

1. 24 

2. 00 

2. 20 

3.55 

695 

360 

. 3970 






689 

120 

.1218 







240 

. 2752 

1800 

1. 145 

1.85 

2. 03 

3. 28 

695 

360 

. 3970 






686 

60 

. 0502 

| 






300 

. 3468 

: 1750 

1.156 

1. 867 

2. 05 

3. 31 

696 

300 

. 3000 


1. 00 

1. 617 

1. 77 

2. 86 

696 

300 

. 3000 






688 

240 

. 2500 







60 

. 0500 

1755 

.83 

1. 34 

1.47 

2. 37 

696 

300 

. 3000 






685 

180 

.1853 







120 

.1147 


. 875 

1.412 

1.55 

2. 50 

696 

300 

. 3000 






687 

175 

. 1676 







125 

. 1324 

1705 

1.06 

1. 71 

1. 87 

3. 02 
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Specimen Dia, 
Specimen Mat'l 
Nozzle Exit Dia. 
Stag'n Press. , atm 
Enthalpy Ratio 
Reynolds No., Re 


TABLE 17 

SUMMARY OF MASS RATE DATA 
1/4” Notes - 

Mat'l. X ( 1 ) Unmodified data 

5” _3 (2) Data modified for 

8. 25 x 10 Reynolds No. Effects 

350 
11 


0 Symbol on graphs 


Run 

Nos. 


Time 

Sec. 


Mass Loss 
Lb/ ft. 2 


Avg. 

°K 

Temp. 


Mass Rate x 10' 
lb/ft 2 - sec 


(1) 


ITT 


‘VI 


Rb 

Pe 


" x 10° 


lb/ft 3/2 -sec-Atm 1/2 


( 1 ) 


(2) 


696 

689 


696 

686 


300 
1 20 
180 

300 

60 

240 


, 3000 
, 1218 
1782 

3000 

0502 

2498 


1775 


1725 





Specimen Dia. 
Specimen Mat'l 
Nozzle Exit Dia. 
Stag'n Press. , atm 
Enthalpy Ratio 
Reynolds No. , Re 


TABLE 18 

SUMMARY OF MASS RATE DATA 

1 / 4 n Notes - 

Mat’l. X (1) 

1 . 2 " ( 2 ) 

5. 60x10' 

350 O 

85 


Unmodified data 
Data modified for 
Reynolds No. Effects 

Symbol on graphs 













Specimen Dia. 
Specimen Mat'l 
Nozzle Exit Dia, 
Stag’n Press. , atm 
Enthalpy Ratio 
Reynolds No. Re 






TABLE 18 

SUMMARY OF MASS RATE DATA 
1/4” Notes - 


Mat'l. X 

1 . 2 ” 

5. 6u x 10 " 2 

350 

85 


(1) Unmodified data 

(2) Data modified for 
Reynolds No. Effects 

© Symbol on graphs 


Run 

Nos. 

Time 

Sec. 

Mass Loss 
Lb/ ft. 2 

Avg. 

° K 

3 

Mass Rate 3c 10 
lb/ft 2 - sec 

• XS 3 

m \[ Pe x 

lb/ft 3/2 -sec-Atm 1/2 

Temp. 

(1) 


(1) 

(2) 

725 

280 

. 3216 






724 

160 

. 1637 







120 

. 1579 


13. 17 

12. 33 

8. 92 

8. 37 

725 

280 

. 3216 






723 

121 

. 1013 







159 

. 2203 

2080 

13.87 

12. 98 

9. 30 

8. 70 

725 

280 

. 3216 






722 

90 

. 0749 







190 

. 2467 

2095 

13. 00 

12. 08 

8. 71 

8. 09 

725 

280 

. 3216 






721 

40 

. 0438 







240 

. 2778 

2090 

11. 58 

10. 76 

7. 84 

7. 30 

725 

280 

. 3216 






728 

40 

. 0515 







240 

1 

. 2701 


11. 26 

10. 77 

7. 62 

7. 29 

726 

240 

. 2707 

1990 

11. 28 

10. 98 

7. 64 

7.45 

726 

240 

. 2707 






724 

1 60 

.1637 







80 

.1070 

2010 

13. 38 

12.92 

9. 06 

8. 76 

726 

240 

. 2707 






723 

121 

.1013 







119 

.1694 

2010 

14. 26 

13. 78 

9.65 

9. 33 

726 

240 

. 2707 






722 

90 

. 0749 







150 

. 1958 

2025 

13. 04 

12. 53 

8.84 

8. 50 

726 

240 

. 2707 






721 

40 

. 0438 







200 

. 2269 

2025 

11.45 

11.00 

7. 75 

7.45 


^SSIFIU 
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TABLE 18 

SUMMARY OF MASS RATE DATA 


Specimen Dia. 
Specimen Mat'l 
Nozzle Exit Dia. 


1/4” 

Mat* 1. X 
1 . 2 " 


Stag'n Press. , atm 5. 60 x 10“ 2 
Enthalpy Ratio 350 

Reynolds No. Re 85 


Notes - 

( 1 ) Unmodified data 

(2) Data modified for 
Reynolds No. Effects 

G Symbol on graphs 


Run Time 
Nos. Sec. 


Mass Loss 
Lb/ ft. 2 


Avg. 

°K 

Temp. 


Mass Rate k 10' 
lb/ft 2 - see 

hi — i — (ZT~ 


m \| Pe * 10* 

lb/ft^-sec-Atm 1 







• • • • 



Specimen Dia. 

1" 

Notes - 


Specimen Mat'l 

Mat' 1. X 

(1) 

Unmodified data 

Nozzle Exit Dia. 

5" 

(2) 

Data modified for 

Stag‘n Press. , atm 

8. 25 x 10' J 


Reynolds No. Effects 

Enthalpy Ratio 

350 

A 

Symbol on graphs 

Reynolds No. , Re^ 

44 



Run 

Nos. 

Time 

Sec. 

Mass Loss 
Lb/ ft. 2 

Avg. 

°K 

Temp. 

Mass Rate k 10^ 
lb/ft^ - sec 

• .1 ?k 3 

m \! Pe x 10 J 

lb/ft 3/2 -sec-Atm 1/2 

a) 

T2J 

(1) 

(2) 

682 

300 

. 1562 






678 

240 

. 1193 







60 

. 0369 

1610 

.615 

. 733 

2. 17 

2.48 

678 

240 

. 1193 

1580 

.499 

. 594 

1.76 

2. 095 

678 

240 

. 1193 






681 

120 

. 0454 







120 

. 0739 

1560 

. 615 

.731 

2. 17 

2. 58 

678 

240 

. 1193 






667 

180 

. 0817 







60 

. 0376 

1565 

. 627 

.746 

2. 21 

2. 63 

678 

240 

. 1193 






679 

180 

. 0763 







60 ! 

1 

.0430 j 

1640 

. 717 

. 853 

2. 52 


679 

180 

. 0763 


.425 

. 506 

1. 50 

1.785 

679 

180 

. 0763 






681 

120 

. 0454 







60 

. 309 

1555 

.602 

. 717 

2. 12 

2. 52 

680 

360 

. 1964 


. 546 

. 650 

1.92 

2. 29 

680 

360 

. 1964 






679 

180 

. 0763 




i 



180 

. . 1199 

1555 

.665 

. 791 

2. 34 

2. 79 

680 

360 

. 1964 






678 

240 

.1193 







120 

. 0771 

1595 

. 642 

. 765 

2. 26 

2. 69 

680 

360 

. 1964 






683 

60 

. 0279 







300 

.1685 

1555 

. 562 

. 670 

1. 97 

2. 34 
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TABLE 20 

SUMMARY OF MASS RATE DATA 


Specimen Dia. 

1 /1 6” 

Notes - 

Specimen Mat'l 

Mat' 1. X 

(1) 

Nozzle Exit Dia. 

5" 

(2) 

Stag'n Press. , atm 

7. 70 x 10’ 3 


Enthalpy Ratio 

350 

• 

Reynolds No. , Re g 

2. 7 



Unmodified data 
Data modified for 
Reynolds No. Effects 

Symbol on graphs 


Run 

Nos. 

Time 

Sec. 

Mass Loss 
Lb/ft. 2 

Avg. 

°K 

Temp. 

3 

Mass Rate k 10 
lb/ft^ - sec 

• S™ 3 

m \| p e x 10 

lb/ ft^-sec-Atm*^ 

a) 

( 2 ) 

(i) 

( 2 ) 

676 

42 

. 0310 

1580 

0. 720 

1. 635 

. 619 

1.405 

674 

120 

. 1594 






675 

60 

. 0758 



i 

1 




60 

. 0836 

1715 

1. 39 

3. 06 

1. 27 

2 . 88 

670 

300 

. 3520 






672 

60 

. 0435 







240 

. 3085 

1690 

1. 29 

2. 93 

1 . 18 

2 . 68 

673 

240 

. 2365 






677 

120 

. 1260 







120 

. 1105 

1670 

0 . 920 

2 . 09 

. 842 

1. 92 

673 

■'.40 

. 2365 






671 

1 30 

. 1661 







60 

. 0704 

1705 

1.17 

2 . 66 

1.06 

2.41 

671 

180 

.1661 

1750 

0 . 92 

2 . 09 

. 842 

1.915 

671 

180 

.1661 






676 

42 

. 0310 







1 38 

. 1351 

1750 

0 . 98 

2 . 22 

. 896 

2. 035 

671 

180 

. I66l 






672 

60 

. 0435 







120 

.1226 

1675 

1.02 

2 . 315 

.933 

2 . 12 

671 

180 

.1661 






677 

120 

.1260 







60 

. 0401 

1715 

670 

1. 52 

.. 613 

1. 39 

677 

120 

.1260 






676 

42 

. 0310 







78 

. 0950 

1630 

1. 22 

2. 77 

1.12 

2. 54 
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Specimen Dia. 
Specimen Mat'l 
Nozzle Exit Dia* 
Stag'n Press. , atm 
Enthalpy Ratio 
Reynolds No. , Re 



TABLE 20 

SUMMARY OF MASS RATE DATA 
1/1 6 M Notes - 

Mat'l. X (1) 

5 " ( 2 ) 

7. 70 x 10 

350 • 

2. 7 


Unmodified data 
Data modified for 
Reynolds No. Effects 

Symbol on graphs 




. 0435 

. 0435 
. 0310 
. 0125 

. 3520 

. 3520 
. 0310 
. 3210 

. 3520 
. 1260 
. 2260 

. 3520 
. 1 661 
. 1859 

. 3520 
. 2365 
.1155 

. 2365 

. 2365 
. 0310 
. 2055 

. 2365 
. 0435 
. 1930 
















TABLE 21 

SUMMARY OF MASS RATE DATA 



1 M 

Mat'l. X 
5 ” 

7.00 x 10 
160 
51 


Notes - 

(1) Unmodified data 

(2) Data modified for 
Reynolds No. Effects 

^ Symbol on graphs 







TABLE 22 MONOTONICALLY SHRINKING SPHERES 
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FIGURE 1. OXIDATION OF GRAPHITE IN NEAR FREE MOLECULE 

HYPERSONIC COUETTEFLOW 
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FIGURE 4. HYPERSONIC ARC HEATED WIND TUNNEL 



CHARACTERISTICS OF HYPERSONIC ARC WIND TUNNEL 








202 


UNCLASSIFIED 




PLASMA JET 



UNCLASSIFIED 


DUMP \ 

GASES 1 'THROTTLING VALVE 









• • 
• • 


UN£t:AS’StFf E'D 



FIGURE 10. TYPICAL MASS TRANSFER MODEL, CALORIMETER 



FIGURE 11. TYPICAL SPHERICAL MODEL CALORIMETERS 
(BEFORE AND AFTER EXPOSURE) 
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MODIFICATIONS OF AT.T GRAPHITE SPECIMENS 
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FIGURE 16. MASS RATE SPECIMEN, 1/16" DIAMETER 



FIGURE 17. INTEGRATED MASS 
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1. 0" DIAMETER SPHERE 
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FLOW REGIME CORRELATION FUNCTION 
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FIGURE 23. 




FIGURE 24. FIXED SPHERES, 1.0" DIAMETER - ATJ ABLATED SPECIMENS 
(Support Axis Co-Axial with Flow Direction During Test: Downward 
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FIGURE 25. ROTATED SPHERES, 1. 0" DIAMETER - AT J SHOWN 

(Axis of Rotation Perpendicular to the Flow Direction During Test: Horizontal) 
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FIGURE 68. COMPARISON OF SURFACE APPEARANCE OF MATERIAL "X" AND AT J GRAPHITE 

AFTER EXPOSURE (Flow Downward) 
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FIGURE 73. MONOTONICALLY SHRINKING SPHERE 
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FIGURE 74. ONE DIMENSIONAL CYLINDER 
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FIGURE 75. TWO-DIMENSIONAL SPHERE(Rotatton About Axis 90°to Flow Direction) 
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FIGURE 76. TWO-DIMENSIONAL (Rotation Abont Axis Aligned to Flow Direction) 

267 


UNCLASSIFIED 



foMAib LOSS ~ C^/^ e */oo) 









JJfLCiASSIFJED. 


• • • • • • 

• • • ••• 



UNCLASSIFIED 















UNCLASSIFIED 

















JJKQLfiS.$IFJED. 


♦ • * 

• • • • 


• • 
• •• 


• • • 
• • 


• • 

• • • 



273 


FIGURE 82. 



UNCLASSIFIED 



















UNCLASSIFIED 







.UN£L ASSIF.I EH 



275 


FIGURE 84. 


UNCLASSIFIED 

















FIGURE 85 



UNCLASSIFIED 

?• 


276 



liN&ASSIFIEQ 






UNCLASSIFIED 




UNJCLASSIF.IED 



UNCLASSIFIED 


{ 



FIGURE 88. 


UNCLASSIFIED 











UNCLASSIFIED 












MCLAS.SIFJEQ. 


•U IlUL Cl %M%J I I h L hP* " • •• 

• •• • • • ••• • •• •• •• 

••••• • ••• *•••••••• 

• • • • • ••• • • • • •• • • 



UNCLASSIFIED 




41 N-G LA$ Sill ED •“ 
.*• *••• •• 



UNCLASSIFIED 






lipASSra.-:; 
•• •• 



FIGURE 92. 


UNCLASSIFIED 









•Up-ASSkNEJ-:; 
• !*• •• •• 



FIGURE 93. 


UNCLASSIFIED 




























; -U NC L A § S tFJ EO • 

•• •• 



UNCLASSIFIED 















• • 


nr-Wl^SIFIED 




290 

UNCLASSIFIED 


UNCLASSIFIED. 


•• 

► • • 


• • * 
• •• • • 



UNCLASSIFIED 












UNCLASSIFIED 

u r | \ j 

• • 



UNCLASSIFIED 







FIGURE 102 



293 


UNCLASSIFIED 







UNCLASSIFIED .. 

: • | 

• • • • • •• •• 



^r:rnrp:pgr 
h:j iftf'-nH' !.ii! 




FIGURE 103. 


UNCLASSIFIED 




















UNCIA.SSLFIEJ).- : 

-r : • : ?. • s " 


»• • • « 
• • 



FIGURE 104. 


UNCLASSIFIED 

















UNCLASSIFIED 








UNCLASSIFIED ... 

• ••• .*• I • • • •, ; 

• * 


• ••• 
• • 

• •• 

• • 

• ••• 


•• 

• • 
• • 



UNCLASSIFIED 























UNCLASSIFIED 

: • 1 ••• > • r L 

• • •• 



UNCLASSIFIED 






UNCLASSIFIED. 


• • • 
• • 

• • 



FIGURE 108. 


UNCLASSIFIED 










UNCLASSIFIED .. 
• ••• /. ! • • • • 

9 # ••• 



FIGURE 109. 


UNCLASSIFIED 





















UNCLASSIFIED 

. : : : % • •• 

* :r :••• ••• •• ' • 


• t • •• 

• • • 

M • * 

• • • 
• •• •• 



FIGURE 110. 


UNCLASSIFIED 















.W.LA.S.SJ-FI.EA 

• • • • • • • J 

• • ••• • • • 




UNCLASSIFIED 












r JUWAjUP : u 

• • •• • ; * 

••• 



UNCLASSIFIED 





.. ...UN CL A.S.S.I-F LE-D: 

• •• • • • • •. \ ! 

• • ••• • • • 

# ••• •• ••• • • •• •• • 



FIGURE 113, 


UNCLASSIFIED 

















WUOTEU . 

^ • ••• • ••• • 
« •• ••• •••••*•• # 



UNCLASSIFIED 







UNCLASSIFIED 


> • • • • • 



FIGURE 115* 


UNCLASSIFIED 














• • 

• • 
• • 




CURE 117 

^CLASSIFIED 


308 















